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ABSTRACT 


Polyphase continuous waveform (CW) radar systems often use the popular Frank code 
and P4 code due to their linear time-frequency characteristics as well as their low 
periodic ambiguity sidelobes. The phase relationship of the Frank code corresponds to a 
sawtooth folding waveform. The phase relationship of the P4 code is symmetrical with a 
parabolic distribution. The radar system’s unambiguous target detection range is limited 
by the number of subcodes within the code period (code length). Increasing the code 
length to extend the unambiguous range results in a larger range-Doppler correlation 
matrix processor in the receiver, a longer compression time and an increase in the 
receiver’s bulk memory requirements. In addition, the entire code period may not be 
returned from the target due to a limited time-on-target resulting in significant correlation 
loss. To significantly extend the unambiguous range beyond a single code period, this 
thesis explores the relationship between the polyphase codes (Frank and P4) 
and the number theoretic transforms (NTT) where the residues exhibit the same 
distribution as the polyphase values. The unambiguous range is extended from the 
number of subcodes within a single code period to the dynamic range of the transform 
without requiring a large increase in correlation processing. The dynamic range of a NTT 
is defined as the greatest length of combined phase sequences that contain no ambiguities 
or repeated paired terms. By transmitting N>2 coprime code periods, the unambiguous 
range can be extended by considering the paired values from each sequence. A 
new Frank phase code formulation is derived as a function of the residue number system 
(RNS) where each residue corresponds to a phase value within the code period (modulus) 
sequence. Based on the symmetrical distribution of the P4 code, a new phase code 
expression is derived using both the symmetrical number system (SNS) and the robust 
symmetrical number system (RSNS). Here each phase value within the code period 
corresponds to a symmetrical residue. MATLAB simulations are used to verify the new 
expressions for the RNS, SNS and RSNS phase codes. Implementation considerations of 
the new approach are also addressed. 


V 



THIS PAGE INTENTIONALLY LEET BLANK 


VI 



TABLE OF CONTENTS 


I. INTRODUCTION.I 

A. CONTINUOUS WAVE RADAR SYSTEMS USING POLYPHASE 

MODULATION.I 

B. PRINCIPAL CONTRIBUTIONS.5 

C. THESIS OUTLINE.6 

II. PRINCIPLES OF PHASE MODULATION RADAR.7 

A. CW PHASE MODULATION CONCEPTS.7 

B. CODE COMPRESSION OF RECEIVED SIGNAL.9 

C. PERIODIC AUTOCORRELATION FUNCTION.14 

D. PERIODIC AMBIGUITY FUNCTION.16 

HI. FRANK POLYPHASE MODULATION AND THE RESIDUE NUMBER 

SYSTEM.19 

A. FRANK PHASE CODE.19 

B. THE RESIDUE NUMBER SYSTEM.22 

C. RESIDUE-FRANK PHASE CODE.24 

D. RESOLVING RANGE AMBIGUITIES USING N RESIDUE-FRANK 

PHASE CODE SEQUENCES.26 

1. Block Diagram of the Residue-Frank Radar System.26 

2. Transmitted and Reference Codes for Compression.27 

3. Calculating the Target Range.31 

4. Resolving Multiple Target Range Ambiguities.34 

5. Practical Considerations.36 

a. Detection of Target in Noise Using RNS Compression . 36 

b. Required CW Signal Power Versus Required Output SNR ..41 

c. Unambiguous Detection Range for a Constant Output 

SNR . 43 

IV. P4 POLYPHASE MODULATION AND THE SYMMETRICAL NUMBER 

SYSTEM.45 

A. P4 PHASE CODE.45 

B. THE SYMMETRICAL NUMBER SYSTEM.47 

C. SYMMETRICAL RESIDUE-P4 PHASE CODE.50 

D. RESOLVING RANGE AMBIGUITIES USING N SYMMETRICAL 

RESIDUE-P4 PHASE CODE SEQUENCES.51 

1. Block Diagram of the Symmetrical Residue-P4 Radar System.51 

2. Transmitted and Reference Codes for Compression.52 

3. Calculating the Target Range.55 

4. Resolving Multiple Target Range Ambiguities.59 

5. Practical Considerations.60 

a. Detection of Target in Noise Using SNS Compression . 60 

b. Required CW Signal Power Versus Required Output SNR ..65 

vii 



































c. Unambiguous Detection Range for a Constant Output 

SNR . 66 

V. P4 POLYPHASE MODULATION AND THE ROBUST SYMMETRICAL 

NUMBER SYSTEM.69 

A. P4 PHASE CODE.69 

B. THE ROBUST SYMMETRICAL NUMBER SYSTEM.69 

C. ROBUST SYMMETRICAL RESIDUE-P4 PHASE CODE.72 

1. Example of Phase Sequence Calculation.73 

2. ACE, PACE, and PAF.74 

D. RESOLVING RANGE AMBIGUITIES USING N ROBUST 

SYMMETRICAL RESIDUE-P4 PHASE CODE SEQUENCES.79 

1. Block Diagram of the Robust Symmetrical Residue-P4 Radar 

System.79 

2. Transmitted and Reference Codes for Compression.79 

3. Calculating the Target Range.83 

4. Resolving Multiple Target Range Ambiguities.86 

5. Practical Considerations.87 

a. Detection of Target in Noise Using RSNS Compression . 87 

b. Plot of Required Signal Power Versus Required SNR 

Output . 92 

c. Unambiguous Detection Range for a Constant Output 

SNR . 93 

VI. CONCLUDING RESULTS AND RECOMMENDATIONS.97 

A. EQUATION SUMMARY.97 

B. COMPARISON OF 3-CHANNEL MODULAR NUMBER SYSTEMS 

WITH THE SAME DYNAMIC RANGE.99 

C. COMPARISON OF 3-CHANNEL MODULAR NUMBER SYSTEMS 

WITH THE SAME MODULI.100 

D. RANGE DETECTION ERROR CONTROL.lOI 

E. RECOMMENDATIONS FOR FUTURE WORK.lOI 

APPENDIX - LIST OF VARIABLES.103 

LIST OF REFERENCES.105 

INITIAL DISTRIBUTION LIST.109 


viii 




























LIST OF FIGURES 


Figure 1. Geometry of CW radar using phase modulation eodes to sample a moving 

target. 8 

Figure 2. Phase eoded signal with the Barker eode N^=7, A = l, /,=! kHz , 

t^= \ ms.9 

Figure 3. The relationship between uncompressed pulse train and phase-coded CW 
signal for (a) Transmitted signal, and (b) The output of the signal 

processed in the receiver. 10 

Figure 4. Correlation receiver matched to periods of a transmitted polyphase 

code for V = 0 (After [9]).12 

Figure 5. Example of code compression in CW Barker code N^=l and N ^ =3.13 

Figure 6 . Normalized correlation output of Barker code sequence N^=7 and 

= 1 for (a) ACF and (b) PACF.15 

Figure 7. PAF of Barker code N^=7 with (a) N^= \ and (b) =4 .17 

Figure 8 . Frank code modulation for M = 8 , = 64 (From [1]).20 

Figure 9. Signal phase (radians) modulo In versus k -index for phase change of the 

Frank code with M = 8 , = 64 (From [1]).20 

Figure 10. ACF (a) and PACF (b) for the Frank code with M = 8 , cpp = 1, and A^ = 1 

(From [1]).21 

Figure 11. PAF for the Frank code modulation with M = 8 , cpp = 1, and A^ = 1 

(From [1]).21 

Figure 12. The RNS residue folding waveforms for m, = 4 and mj = 5.23 

Figure 13. Block diagram for the radar using A residue-Frank phase code sequences...26 

Figure 14. Illustration of transmitted signal using A = 3 residue-Frank phase code 

sequences for m, = 7>,m^= 4, and =5.28 

Figure 15. Plot of the residue-Frank phase code sequence for mj = 3, = 4, and 

m 3 =5 (M ^^5 =60) showing (a) residues, and (b) phase sequences from 

(22).29 

Figure 16. Illustration of compression at the receiver and the range bin for RNS 

m, = 3, A^ = 9, and A^ = 1.30 

Figure 17. Example of compression output for target detection using A = 3 residue- 

Erank phase code sequences with m, =1), m^= 4, and m^=5 .33 

Eigure 18. Resolving the range ambiguity of 2 targets using RNS for 

m, = 3, m 2 = 4, and m 3 = 5.35 

Eigure 19. Missed detections of the first set of residues.35 

Eigure 20. Residue-Erank signal without noise for m, = 3, m^= 4, and m^=5 with 

^RNs = 60, A = 1, /^ = 1 MHz, /^ = 7 MHz, cpp = 1, and C = 1 jUS .36 






















Figure 21. Residue-Frank signal with SNR = 30 dB for m, =3, =4, and m^=5 

with = 60, A = 1, /^ = 1 MHz, /, = 7 MHz, cpp = 1, and = 1 ps .37 

Figure 22. Residue-Frank signal with SNR = 0 dB for m, = 3, m 2 = 4, and m^=5 with 

Mj^Ns = 60, A = 1, /^ = 1 MHz, f^-1 MHz, cpp = 1, and = 1 //s.37 

Figure 23. Target deteetion for SNR = 30 dB showing the range bins for 

m, = 3, m 2 = 4, and m^=5 with the two targets at 440 m and 5840 m.39 

Figure 24. Target deteetion for SNR = 0 dB showing the range bins for 

mj =3, 1712= 4? and m 3 = 5 with two targets at 440 m and 5840 m.40 

Figure 25. Illustration of deteetion error from the threshold deteetor (CFAR) using 

the residue-Frank for m, = 3, m 2 = 4, and =5.41 

Figure 26. Average power of the CW transmitter for residue-Frank with 

m, = 3, m 2 = 4, and m 3 = 5.42 

Figure 27. Comparison of the maximum unambiguous range of CW radar system for 
SNRf,^ = 13 dB using the residue-Frank phase eode for = 3, m 2 = 4, and 
m 3 = 5 with the radar using eaeh Frank eode sequenee individually with 

the eorresponding .44 

Figure 28. P4 phase sequenee for = 64 (From [1]).46 

Figure 29. Signal phase (radians) modulo Itt versus k -index for phase ehange of the 

P4 eode with = 64 (From [1]).46 

Figure 30. ACF (a) and PACF (b) for the P4 eode with = 64, cpp = 1,and N^=\ 

(From [1]).47 

Figure 31. PAF for the P4 eode modulation with = 64, cpp = 1, and N^=\ 

(From [1]).47 

Figure 32. The SNS residues for m, = 4 and m 2 =5 .48 

Figure 33. Bloek diagram for the radar using N symmetrieal residue-P4 phase eode 

sequenees.52 

Figure 34. Illustration of transmitted signal using N = 3 symmetrieal residue-P4 

phase eode sequenees for = 7, m 2 = 8 , and = 9.53 

Figure 35. Plot of the symmetrieal residue-P4 for mj = 7, m 2 = 8 , and m 3 =9 

(^sNs =37) showing (a) symmetrieal residues, and (b) phase sequenees 

from (34).54 

Figure 36. Illustration of eompression at the reeeiver and the range bin for SNS 

m, = 7, = 7,and N^=\ .55 

Figure 37. Illustration of target deteetion by using the range bin matrix for SNS-P4 

eode with = 7, m 2 = 8 , and = 9.57 

Figure 38. Resolving the range ambiguity of 2 targets using SNS for 

m, = 7, m 2 = 8 , and = 9 .59 

Figure 39. Missed deteetions of the first set of symmetrieal residues.60 


X 






















Figure 40. Symmetrical residue-P4 signal without noise for m, = l,m^= 8 , and 
m 3 = 9 with = 37, A = 1, /^ = 1 MHz, f^=l MHz, cpp = 1,and 

t^= \ ps .61 

Figure 41. Symmetrical residue-P4 signal with SNR = 30 dB for m, = 7, = 8 , and 

m 3 = 9 with = 37, A = 1, /^ = 1 MHz, f^=l MHz, cpp = 1, and 

t^= \ ps .61 

Figure 42. Symmetrical residue-P4 signal with SNR = 0 dB for m, = l,m^= 8 , and 

m 3 = 9 with = 37, A = 1, /^ = 1 MHz, = 7 MHz, cpp = 1, and 

= 1 ps .62 

Figure 43. Target detection for SNR = 30 dB showing the range bins for 

m, = l,m^= 8 , and m^=9 with two targets at 890 m and 2990 m.63 

Figure 44. Target detection for SNR = 0 dB showing the range bins for 

m, = 7, m 2 = 8 , and m^=9 with two targets at 890 m and 2990 m.64 

Figure 45. Illustration of detection error from the threshold detector (CFAR) using 

the symmetrical residue-P4 for = l,m^= 8 , and = 9.65 

Figure 46. Average power of the CW transmitter for symmetrical residue-P4 with 

m, = 7, m 2 = 8 , and = 9. 66 

Figure 47. Comparison of the maximum unambiguous range of a CW radar system 
for SNRj^^ =13 dB using symmetrical residue-P4 phase code for =7, 

= 8 , and m^=9 with the radar using each P4 code individually with the 

corresponding .67 

Figure 48. RSNS residues RS^, for modulus m = 16, A = 2, and =64.75 

Figure 49. The robust symmetrical residue-P4 phase code sequence using m = 16, 

N = 2, = 64 is shown with the relationship with the P4 code using 

=32.76 

Figure 50. Power spectral density of the robust symmetrical residue-P4 signal using 

m = 16, A = 2, and A^ = 64 with f, = 1 kHz, / = 7 kHz, and cpp = 2 .77 

‘^RSNS ^ 

Figure 51. The ACF and PACT of RSNS phase-coded signal 
(m = 16, A = 2, A^ =64) at the receiver with A =1, 

^ ^RSNS ^ ^ 

/, = 1 kHz, /j = 7 kHz, and cpp = 2.77 

Figure 52. The PAF of the signal for (a) the robust symmetrical residue-P4 phase 
code sequence using m = 16, A = 2, = 64, and tl= \ ms, and (b) the 

P4 code using A^ = 32, and C = 2 ms.78 

Figure 53. Block diagram for the radar using A robust symmetrical residue-P4 phase 

code sequences.79 

Figure 54. Illustration of transmitted signal using A = 3 robust symmetrical residue- 

P4 phase code sequences for m, =2,m^= 4, and = 5.81 


XI 


















Figure 55. Plot of the RSNS using m, =3, =4, and =5 with =43 for 

(a) folding waveforms, and (b) their phase sequences.82 

Figure 56. Illustration of compression at the receiver and the range bin for RSNS 

m, = 3, = 18, and = 1.83 

Figure 57. Illustration of target detection by using the range bin matrix for RSNS-P4 

code with =3, m^= 4, and m^=5 .85 

Figure 58. Resolving the range ambiguity of 2 targets using RSNS for 

m, = 3, m 2 = 4, and m 3 = 5. 86 

Figure 59. Robust symmetrical residue-P4 signal without noise for 

m, = 3, m 2 = 4, and m^=5 with 

^RSNS = 43, A = 1, /^ = 1 MHz, = 9 MHz, cpp = 1, and tl = 0.33 ps .87 

Figure 60. Robust symmetrical residue-P4 signal with SNR = 30 dB for 

m, =3, m 2 = 4, and m 3 = 5 with 

^RSNS = 43, A = 1, /, = 1 MHz, = 9 MHz, cpp = 1, and tl = 0.33 ps . 88 

Figure 61. Robust symmetrical residue-P4 signal with SNR = 0 dB for 

m, = 3, m 2 = 4, and m 3 = 5 with 

Mrsns = 43, a = 1, /, = 1 MHz, = 9 MHz, cpp = 1, and tl = 0.33 ps . 88 

Figure 62. Target detection for SNR = 30 dB showing the range bins for 

m, = 3, m 2 = 4, and m 3 = 5 with two targets at 140 m and 1040 m.90 

Figure 63. Target detection for SNR = 0 dB showing the range bins for 

m, = 3, m 2 = 4, and m 3 = 5 with two targets at 140 m and 1040 m.91 

Figure 64. Illustration of detection error from the threshold detector (CFAR) using 

the robust symmetrical residue-P4 for m, = 3, m 2 = 4, and m 3 = 5.92 

Figure 65. Average power of the CW transmitter for robust symmetrical residue-P4 

with m, = 3, m 2 = 4, and = 5.93 

Figure 66 . Comparison of the maximum unambiguous range of CW radar system for 
SNRf,^ =\3dB using robust symmetrical residue-P4 code for m^=3, 
m 2 = 4, and m 3 = 5 with the radar using each P4 code individually for the 
corresponding N=N^ / N .94 

^ ^ c CrSNS 


xii 















LIST OF TABLES 


Table 1. Finding the RNS dynamic range for m^=A and m 2 =5 .23 

Table 2. Comparison of the residue-Frank phase code sequence (m = 4) with the 

Frank phase code sequence {M =4, =\6) .25 

Tables. Bin matrix and range intervals of 3-channel RNS with 

m, = 3, m 2 = 4, and m 3 = 5.31 

Table 4. Finding the SNS dynamic range for m, = 4 and mj = 5.49 

Table 5. Comparison of symmetrical residue-P4 phase code sequence 

(m = 8, = 8) with the P4 phase code sequence (N^=8) .51 

Table 6 . Bin matrix and range intervals of 3-channelSNS for m, = 7, mj = 8 ,and 

m 3 = 9.56 

Table 7. The RSNS folding waveforms for m, = 3 (i'j = 0), mj = 4 (S 2 = 1), and 

m 3 = 5 ( 5'3 = 2) (After [20]).72 

Table 8 . Comparison of the robust symmetrical residue-P4 phase code sequence 

(m = 4, N = 2, =16) with the P4 phase code (N=8) .73 

Table 9. Bin matrix and range intervals of 3-channel RSNS for 

m, = 3, m 2 = 4, and = 5.84 

Table 10. Summary of equations for each modular phase code.97 

Table 11. Summary of equations for N - channel modular number system.98 

Table 12. Comparison of each modular phase code in which each has the same 

dynamic range = 60, and tl=\ jus .99 

Table 13. Comparison of each modular number system in which each has the same 

moduli [3 4 5], and = ljus .100 

Table 14. List of variables.103 


xiii 

















THIS PAGE INTENTIONALLY LEET BLANK 


XIV 



EXECUTIVE SUMMARY 


The Frank phase code and the P4 phase code are often used in continuous wave (CW) 
polyphase modulation radar systems to detect the target’s range due to their linear time- 
frequency characteristics as well as their low periodic ambiguity sidelobes. Both phase 
codes limit the unambiguous range to within a single code period T = , where 

is the number of subcodes and is the subcode period. To extend the unambiguous 

detection range, the radar might simply increase the code period of the polyphase 
modulation signal. However, the longer code period requires a larger range-Doppler 
correlation matrix processor, an increase in the compression time, as well as an increase 
in the bulk memory in the receiver. Moreover, a significant amount of loss is incurred if 
the total number of code periods returned from the target is less than the number of 
reference code periods used in the compression process due, for example, to a limited 
time-on-target. A new set of waveform architectures that can solve these limitations is 
proposed in this thesis. 

Number theoretic transform s (NTT) based on residue number systems are used in 
many digital signal processing applications to increase the amount of information 
available from various folding waveforms. With a set of A > 2 moduli, the paired integer 
residue terms from the transform sequences must be distinct (no ambiguities) throughout 
the entire dynamic range. The dynamic range of a NTT is defined as the greatest length 
of combined sequences of integers that contain no ambiguities or repeated paired terms. 
Consequently, the NTT residues within the dynamic range can be used to construct a 
low-power CW polyphase waveform that is useful for detecting targets beyond the single 
code period unambiguous range. 

To significantly extend the unambiguous range beyond a single code period, this 
thesis explores the relationship between the polyphase codes (Frank and P4) and the 
NTT where the residues exhibit the same distribution as the polyphase values. The 
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unambiguous range is extended from the number of subcodes within a single code period 
to the dynamic range of the transform without requiring a large increase in correlation 
processing. 

A new Frank phase code formulation is derived as a function of the residue 
number system (RNS) where each residue corresponds to a phase value within the code 
period (modulus) sequence. Based on the symmetrical distribution of the P4 code, a 
new phase code expression is derived using both the symmetrical number system (SNS) 
and the robust symmetrical number system (RSNS). Here each phase value within the 
code period corresponds to a symmetrical residue. MATLAB simulations are used to 
verify the new expressions for the RNS, SNS and RSNS phase codes. Implementation 
considerations of the new approach are also addressed. 

The relationship between the Frank phase code and the RNS is developed based 
on their similar saw-tooth waveform distribution. A new equation using the RNS to 
derive the Frank phase code is formulated. This new phase code sequence is called the 
residue-Frank phase code. The periodic autocorrelation function (PACT) and the 
periodic ambiguity function (PAF) show that the new residue-Frank phase code formula 
gives exactly the same results as using the formula in the Frank phase code. Next, a set of 
N = 3 RNS coprime moduli [3 4 5] (with dynamic range = 60) is chosen to 
examine the idea of extending the unambiguous range. A MATLAB program is used to 
plot the results for detecting two targets that lie in an ambiguous detection range bin if a 
single modulus is used. To solve the unambiguous range, the combined residues from N 
sequences are used. The result shows that the paired terms from both targets are different, 
meaning that both targets have different corresponding ranges. Thus, the residue-Frank 
phase code sequences can be used to extend the unambiguous detection range from 
K =cT/2 = ct,NJ2 to = ct^Mj^j^^ / 2 where c is the speed of light in free space. 

A block diagram of a polyphase CW radar system that uses the residue-Frank phase code 
is shown below in Figure 1. 
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Figure 1. Block diagram of a CW polyphase radar system that uses a residue-Frank 
phase code. 


The relationship between the P4 phase code and the SNS is developed based on 
their similar parabolic waveform distribution. A new equation using the SNS to derive 
the P4 phase code is formulated. This new phase code sequence is called the symmetrical 
residue-P4 phase code. The PACT and the PAF show that the symmetrical residue-P4 
phase code formula gives the same results as using the formula in the P4 phase code. 
Then, the symmetrical residue-P4 phase code sequences using a set of A = 3 SNS 
coprime moduli [7 8 9] (with dynamic range M^ns = 37 ) are examined for the two targets 
located in an ambiguous range bin by using MATLAB programming. Similar to the 
residue-Frank phase code sequences, the result shows that the symmetrical residue-P4 
phase code sequences can be used to extend the unambiguous detection range from 

K=cT 12 = ct,N^ / 2 to = ct,M / 2. 

Reconstruction of the target’s range from the paired RNS and SNS phase terms 

can contain inaccuracies when the incoming target range straddles two range bins. As a 

result, the recovered range of the target has a large error. To prevent the target range 

error, the RSNS is used. The paired term s from all N sequences in the RSNS when 

considered together change one at a time at each range bin subcode transition (integer 

Gray code property). This property makes the RSNS particularly attractive for controlling 
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the target range error. In this thesis, the relationship between the P4 phase eode and the 
RSNS is also developed. A new equation using the RSNS to derive the P4 phase eode is 
formulated. This new phase code sequence is called the robust symmetrical residue-P4 
phase code. Within each P4 subcode, the phase shifts in the robust symmetrical residue- 
P4 phase sequences are repeated N times (equal to the number of moduli), resulting in 
the broader mainlobe for the PACF and the reduction of signal bandwidth as a function of 
N. The PACF and the PAF show that the robust symmetrical residue-P4 phase code 
formula gives similar results as using the formula for the P4 phase code. A MATLAB 
program is used to examine the abilities to extend the unambiguous range and to control 
the range bin errors by using a set of N = 3 RSNS coprime moduli [3 4 5] (with dynamic 
range = 43). The results verify that the unambiguous range can be extended from 

R^=cT/2 = ct,NJ2 to = ct[M / 2 where t[ is the subcode period in the 

RSNS phase code, and the range bin has an error of only one range bin. 

In summary, all three NTT are investigated to show that they can be used to 
extend the unambiguous range. Also, the target range error from coding can be prevented 
by using the robust symmetrical residue-P4 phase code sequences. The use of these new 
derived phase sequences can be applied to CW radar systems in order to have a low 
probability of intercept characteristic given that the set of moduli and the radar signal 
parameters are chosen properly. The next step in this research is to investigate the 
performance of using various sets of coprime moduli, and the hardware implementation 
and testing of an actual system and is being proposed for future work. 
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I. INTRODUCTION 


A. CONTINUOUS WAVE RADAR SYSTEMS USING POLYPHASE 

MODULATION 

Continuous wave (CW) radar systems have a superior low probability of intereept 
(LPI) performance over pulse train radar systems due to their low average power 
transmitted (e.g., 1 W) and their use of pulse compression techniques. Note we extend the 
concept of pulse compression to unpulsed CW waveforms since the techniques are 
similar and the objectives are the same. Power management is also a benefit when a solid 
state phased array antenna is used. This enables the target return signal-to-noise ratio 
(SNR) to remain at a constant value within the radar receiver, thus lowering the 
transmitted power further as the range-to-target decreases. An LPI radar system is 
defined as a radar system that uses a special emitted waveform intended to detect targets 
as well as to prevent a non-cooperative intercept receiver from intercepting and detecting 
its emission [1]. On the other hand, pulsed radar systems require a relatively high peak 
power (e.g., 10 kW) to be transmitted to obtain the same probability of target detection. 
Consequently, the non-cooperative intercept range is significantly longer for the pulsed 
radars. 

Since pure CW waveforms cannot resolve the target’s range, periodic modulation 
techniques are used, such as frequency modulated CW (FMCW), frequency-shift keying 
(FSK), noise modulation, phase-shift keying (PSK), as well as hybrids of these 
techniques. The first step to designing CW radar systems that use periodic modulations 
for compression is to decide on the range resolution that is required. This in turn sets the 
transmitted bandwidth of the waveform for the above techniques (except for FSK where 
the range resolution is dependent on the duration of each frequency). Due to the advances 
in high-speed processing and direct digital synthesis modules [2], the use of PSK 
techniques in CW radar is highly advantageous. CW radars that transmit and receive PSK 
signals can result in LPI radar systems with small range resolution cells and are ideally 
suited for many sensor applications for situational awareness, including multiple-input 
multiple-output (MIMO) configurations. 
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Binary Barker phase coded sequences (1, -1) are one of the most popular PSK 
techniques in the design of pulse compression radar signals due to the aperiodic 
autocorrelation coefficients or time sidelobes being limited to 1 / (at zero Doppler) 

relative to a mainlobe level of one, where is the number of subcodes used. The 
longest code length is N^= 13 (or A^^=169 for a compound code) [3]. These codes are not 
used in CW radar systems for LPI applications since they are very sensitive to Doppler 
shifts and can be easily detected by an intercept receiver that uses frequency doubling 
detection, which consists of multiplying the signal by itself and processing the result 
through an envelope detector. 

Polyphase sequences are also of finite length with subcodes and consist of 

discrete time complex values with constant amplitude but with a variable phase. 
Polyphase coding refers to phase modulation of the CW carrier with a polyphase 
sequence consisting of a number of discrete phase states with each phase corresponding 
to a subcode. The number of subcodes is taken from an alphabet of size > 2, which is 
also the processing gain of the radar excluding any post detection integration. Increasing 
the number of phase values in the sequence allows the construction of longer sequences, 
resulting in a greater processing gain or compression ratio in the receiver (or equivalently 
a larger SNR). The code period is T = where the size of each subcode 4 determines 

the 3 dB bandwidth of the waveform B = 1 /and the range resolution AR = cti^ / 2. The 
unambiguous range is limited by the number of subcodes in the sequence as = N^AR. 

With a CW polyphase waveform, the matched filter in the receiver is a coherent, 
range-Doppler matrix correlation processor that performs a cross correlation between the 
received signal and a reference signal whose envelope is the complex conjugate of 

code periods of the transmitted polyphase signal. The cross correlation output values are 
added together to reduce the ambiguity sidelobe levels. When the polyphase signals are 
returned and have an impressed Doppler shift, the correlation process used in the receiver 
compression operation is not perfect, resulting in a certain amount of correlation loss due 
to the phase shifts across the code period T. 
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Two of the most important PSK codes that are used in CW radar systems include 
the Frank code [4] and the P4 code [5]. The Frank code is a polyphase code that is 
derived from a step approximation to a linear FMCW waveform using M frequency 
steps and M samples per frequency, giving the number of subcodes or processing gain as 
=M^. The phase steps of the Frank code exhibit a saw-tooth folding waveform 
distribution. The peak sidelobe level (PSL) of a single compressed Frank code period is 
P5L = 201ogjo(l/(M;7')). The P4 code consists of discrete phases of a linear chirp 

waveform taken at specific time intervals. It is derived by converting a linear FMCW to 
baseband using a local oscillator on one end of the frequency sweep and sampling the 
inphase and quadrature video at the Nyquist rate. The P4 phase steps exhibit a 
symmetrical parabolic distribution. The PSL of a single P4 compressed code period is 

PSL = 20log,, ^2/(Ny). 

Although polyphase CW modulations can serve as a transmission waveform for a 
LPI radar system, there are several limiting factors. Most significantly, the unambiguous 
detection range of the waveform is limited by the number of subcodes (code length) as 
= N^AR. By increasing the number of subcodes to extend the unambiguous range, a 
larger range-Doppler correlation matrix processor in the radar receiver is needed to 
provide the code compression. To lower the peak side lobes, copies of the phase code 

are used to compress the returned code periods from the target. Consequently, the 

increased code length or processing gain requires a larger code compression time as well 
as an increase in the bulk memory requirements in the receiver. Further, a significant 
amount of correlation loss is incurred if the total number of code periods returned from 
the target is less than the number of code periods used in the correlation 

processor, as would be the case due to a limited time-on-target. 

In this thesis, a novel relationship between the Frank code and the residue number 
system (RNS) is developed in order to significantly extend the unambiguous range 
beyond a single code period T = Nji^ while not having to increase the number of 

subcodes. In addition, due to the symmetrical distribution of the P4 polyphase code, the 
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unambiguous range is extended by developing a phase relationship based on symmetrical 
number systems; both the symmetrical number system (SNS) and the robust symmetrical 
number system (RSNS) are investigated. The phase relationships are constructed using 
these number theoretic transforms to extend the unambiguous range of the radar system 
by associating each phase subcode within the code period T with an integer sequence 
within a coprime modulus from the set . The paired integer terms from 

each polyphase sequence are unambiguous within the transform’s dynamic range. The 
dynamic range is defined as the greatest length of combined phase sequences that contain 
no ambiguities or repeated paired terms. The dynamic range for the RNS is denoted 

^RNs (Frank code). The dynamic range for the SNS is (P4 code). The dynamic 

range for the RSNS is (P4 code). 

For each case, the length of the dynamic range contains multiple code periods. In 
the radar receiver, the integer values within each modulus sequence correspond to the 
phase values from a phase detector [6]. The RNS phase sequence contains multiple Frank 
codes and extends the unambiguous target detection range from R^=cT 12 = ctf^N^ / 2 to 
RNS ~ m.) 12 = ctf^M / 2. The SNS phase sequence contains multiple P4 codes 

and extends the unambiguous target detection range beyond a single code period to 
SNS ~ SNS / F is noted that the SNS also has a direct relationship with the 

discrete Fourier transform (DFT) and has been used to resolve frequency ambiguities in 
undersampling digital receivers [7]. Similarly, the unambiguous range using the RSNS is 

extended to ! 2 where t[ is the subcode period in the RSNS phase 

code. 

Reconstruction of the target’s range from the paired RNS and SNS phase terms 
can contain inaccuracies when the incoming target range straddles two range bins or 
subcodes. For example, if the target is straddling a range bin and a constant false alarm 
rate (CFAR) processor detects the target in the wrong subcode within a modulus 
sequence, the recovered range of the target will have a large error. This is shown to be the 
case for both the RNS and the SNS. Additional signal processing can be used to eliminate 
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these errors but adds additional complexity to the detection process. A more efficient 
method to eliminate the target range bin errors uses the RSNS. The RSNS is also a 
modular scheme; however, the paired terms from all N sequences, when considered 
together, change one at a time at each range bin subcode transition (integer Gray code 
property). Although the dynamic range of the RSNS is less than the SNS 

{Mr^ns < ^sNs )’ the RSNS integer Gray code property makes it particularly attractive for 

controlling the target range error. In summary, all three number theoretic transforms are 
investigated to determine the feasibility of extending the unambiguous range of the CW 
polyphase emitter. 

B. PRINCIPAL CONTRIBUTIONS 

The first step in this thesis was to understand the principles of CW phase 
modulation radar systems. This included the complex envelope representation of the 
transmitted signal and the related modulation parameters such as the code period, the 
subcode period, the number of subcodes within the code period, and the number of cycles 
of the carrier frequency per subcode. Their relationship to the signal bandwidth, the range 
bin size, the unambiguous detection range and the processing gain were also studied. 
Upon reception, the CW signal is compressed, and several code compression techniques 
were studied and modeled. To understand the receiver performance when a particular 
CW phase code is used, the periodic autocorrelation function and the periodic ambiguity 
function were examined. A MATLAB code was written to generate the CW phase coded 
signals and to call the periodic ambiguity function subroutines. 

The next step in the thesis was to study the Frank polyphase code formulation and 
the phase distribution within a code period. Periodic ambiguity analysis was also 
performed to quantify the sidelobe levels within the receiver’s ambiguity space (range 
offset - Doppler offset). Next the RNS was studied as a means to generate the Frank code 
since they both have sawtooth folding waveform representations. By using the RNS to 
create the Frank polyphase modulation, the unambiguous range of the code can be 
extended significantly beyond a single code period. A new expression to generate the 
Frank polyphase code was developed as a function of the coprime RNS moduli chosen. 
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By using two or more Frank code sequences that are relatively prime, the detected target 
residue from each sequence is combined within a system of congruences which can be 
solved in a straightforward manner using the Chinese Remainder Theorem (CRT). The 
solution is the target’s unambiguous range. A MATLAB model was developed to 
simulate the new Frank polyphase modulation in order to verify the expression and 
determine the differences in phase values (if any) that might occur. Next a new CW radar 
signal processing algorithm was developed to transmit and receive the coprime Frank 
sequences to demonstrate the extension of the unambiguous range. The algorithm also 
addressed multiple targets within the field of view. Several examples were given to 
demonstrate the concept. To examine the practical considerations, a simulation of the 
transmit power versus required compression output SNR was performed in MATLAB for 
several values of coprime moduli. 

Next, the P4 polyphase code was examined the relationships to the SNS and the 
RSNS. A new formulation of the P4 polyphase code as a function of the SNS and RSNS 
coprime moduli was developed. Analysis similar to that for the Frank code described 
above was performed to verify the new P4 expression. Finally, a summary of the different 
relationships that were developed is given and a comparison was performed. 

C. THESIS OUTLINE 

In Chapter II, the fundamentals of radar systems that use polyphase modulation 
are reviewed. In Chapter III, the Frank code and RNS are examined to extend the 
unambiguous range. In Chapter IV, the relationships between the P4 and the SNS are 
developed. In Chapter V, the relationships between the P4 and the RSNS are investigated 
to reduce the possibility of a range bin encoding error. Concluding remarks are offered in 
Chapter VI. 
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II. PRINCIPLES OF PHASE MODULATION RADAR 


The idea of CW radar using phase modulation is reviewed in this ehapter. First, 
the fundamentals and use of phase modulation in CW radar systems are explained. Next, 
the eoneepts of signal analysis and the teehnique used in the reeeiver’s signal proeessing 
are presented. Also, the system equations for a phase modulation in CW radar are given 
in order to examine the signal’s properties. These equations will be used to determine the 
new phase modulation eode’s performanee in the following ehapters. 

A. CW PHASE MODULATION CONCEPTS 

In Figure 1, the basie CW radar geometry is illustrated. First, the CW radar 
transmits the phase eoded signal s(t) to the target through the transmitting antenna. 
When the signal propagates to the target at a range R, it refleets baek to the CW radar 
with different eharaeteristies. For example, the signal’s amplitude is less and the 
frequeney is shifted in time in aeeordanee with the relative movement between platforms 
along the line-of-sight {Vcos0 in Figure 1) or Doppler shift v. The return signal x(t) is 
then eaptured by the reeeiving antenna. In the radar system’s reeeiver and digital signal 
proeessor, the eompression proeessing output is used to deteet and extraet the target’s 
eharaeteristies. 

In phase eoded CW radar systems, the phase shifting operation is performed in the 
radar’s transmitter, with the timing information generated from the reeeiver-exeiter. The 
transmitted eomplex signal ean be written as [1] 

= ( 1 ) 

where A is the signal amplitude, is the earrier frequeney, and is the time 

dependent phase modulation eode. The inphase I and quadrature Q representations of 
the eomplex signal from the transmitter ean be represented as 

Iit) = Acosi27rfj + (^^it)) ( 2 ) 


and 


Qit) = AsinilTTfj + ( 0 ). 
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Figure 1. Geometry of CW radar using phase modulation codes to sample a moving 

target. 


For one code period T, the CW signal is shifted in phase every subcode with each 
code period consisting of subcodes. Each subcode with phase has a duration of 
which is the subcode period. For a specific code sequence, the code period is 


T = N/,. 


The transmitted signal can be expressed as 




k=l 


for 0<t <T and zero elsewhere. The complex envelope is 


( 4 ) 

( 5 ) 

( 6 ) 


If cpp is the number of cycles of the carrier frequency per subcode, the bandwidth B of 
the transmitted signal is 

L 1 


B = 


cpp 


(V) 
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To understand the concepts of phase modulation, a binary Barker phase code is 
examined [3]. In Figure 2, an example of one code period with N^=l is shown. Here, 
the carrier frequency is /^ = 1 kHz, the subcode period = 1 ms, and from (7) cpp = 1. 
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Figure 2. Phase coded signal with the Barker code N^=l, A = l, /^=1 kHz , 

= 1 ms. 


The received return signal x{t) from the target can be written as 

x{t) = (8) 

where p is the amplitude attenuation coefficient of the received signal due to target 
scattering and propagation attenuation. The two-way roundtrip delay t for a monostatic 
radar system is t = 2R/c, where c is the speed of light in free space. The Doppler 
frequency v due to the target motion is v = 2V cos{d) / Ji where A is the signal 
wavelength and V is the target velocity. 

B. CODE COMPRESSION OF RECEIVED SIGNAL 

Code compression (CC) is a signal processing technique designed mainly to 
minimize the range resolution cell and increase the SNR of the returned signal. The idea 
of CC is similar to that in signal correlation in pulse train radar. As shown in Figure 3 a, 
the CW radar repeats the phase sequence every code period T. A similar pulse train 
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signal would also have a pulse repetition interval PRI = T. Referring to Figure 3b, if the 
pulse width of the pulse train signal is equal to the width of the subcode period in 

the CW phase coded signal, the output’s mainlobe from both correlator and compressor 
would have the same duration (). These CC peaks represent the targets if the 

peaks are higher than the threshold desired in a constant false alarm rate (CFAR) process. 
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Figure 3. The relationship between uncompressed pulse train and phase-coded CW 
signal for (a) Transmitted signal, and (b) The output of the signal processed in the 

receiver. 


The phase-coded signal is transmitted to detect the target at a range R. When 
reaching the target, the reflected signal returns to the radar. At the radar receiver, the 
phase coded return signal is correlated with the reference phase code. This is achieved by 
convolving the incoming signal’s phase code with a conjugated and time-reversed 
version of the transmitted signal’s phase code. This operation can be done either in the 
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digital or analog domain. Many digital signal processing architectures have been 
proposed for fast correlation using fast Fourier transforms. If no compression loss is 
present, the correlated output will give a very narrow pulse width at the compression 
output for the detected target. The correlation output is then sent to a target detection 
process such as a constant false alarm rate processor. If there is any alteration of the 
return signal’s phase (for example from a Doppler shift or a less than adequate bandwidth 
in the receiver B <\/t, ), a certain degree of correlation loss occurs. 

The received waveform from the target is digitized and correlated using a phase 
code compressor. The range-Doppler correlation matrix can contain a cascade of sets 

of reference coefficients to improve the sidelobe structure. In Figure 4, the range- 
Doppler correlation matrix receiver for zero Doppler offset is illustrated. Here, the 
received signal consists of code periods, each of which has subcodes. The return 

signal is first processed by a filter matched to a rectangular subcode of length 4 , 
followed by a phase detector that sends forward the phase values. The detected output 
signal is then sent through a tapped delay line where each delay D is 4 . As the signal 
progresses through the tapped delay line, it is multiplied by the reference signal. At each 
step, the multiplication for each delay is summed separately for each of the reference 

code periods. To reduce the sidelobe structure, a weighting function C. can also be added 
[ 8 ]. To compute the entire range-Doppler ambiguity function, the correlation outputs 
from each subcode are multiplied by < 7 ^ = for v = A v where k ranges from 0 to 

-1. For V = 2Av, the output of the previous u = A u operation are again multiplied 
by < 7 ^ = and so on for the entire ambiguity space [9]. 
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Conelation Output 


Figure 4. Correlation receiver matched to periods of a transmitted polyphase code 

for V = 0 (After [9]). 

In Figure 5, an illustration of the code compression for CW binary Barker code 
with N^=l is shown. The radar receives the returned Barker sequence with the number 

of code periods N^ = 3. The number of reference code periods used in the receiver is 
N^=\. The receiver starts compressing the first returned code 7^ at r = 0, giving the 
correlation output peak at a normalized delay of r / C = 0 representing the detected 
target. At each step in the compression, the returned signal is shifted one subcode period 
r = t^. When the returned code is shifted by one code period T, the correlation output 

gives another peak value due to the second code being compressed at a normalized 
delay of r/t^ =7. The last peak from the third code compressed is also shown at a 
normalized delay offset of r /= 14. That is, the peak value for the detected target 
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repeats every code period T. Generally, the width of the PACF mainlobe is 2t^. To be 
able to distinguish between two targets close to one another, the pulses have to be 
separated by at least r = period (which is defined as the range resolution cell). 
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Figure 5. Example of code compression in CW Barker code N^=7 and N^ = 3. 


In summary, the unambiguous range of a specific CW phase-coded sequence is 

cT cM t 

/?„= —= ^ (9) 

2 2 

and the range resolution for one transmitted polyphase signal is 


A/? = ^. 
2 


( 10 ) 


Another important parameter in pulse compression is the pulse compression ratio 
(PCR) or time-bandwidth product. This is also referred to as the radar’s processing gain 
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(PG) and is the preferred term in this thesis. In general, the higher the PG, the higher the 
SNR improvement in the signal proeessing step. The radar signal proeessing input SNR is 
SNRjfi and the radar signal processing output SNR is 5AT?^^(see Figure 1). They are 
related to the processing gain of the signal processor by 


PG 


SNR,„ 
SNR,, ■ 


( 11 ) 


For a polyphase code, the PG is equal to the time-bandwidth product or the number of 
subcodes as 


PG=Til/t,) = iNj,)/t,=N^. (12) 

In the Barker code N^=l example, the processing gain is 7. For a carrier 
frequency /^, narrowing the subcode period enables the transmitted signal to be 
spread over a larger bandwidth. Also, increasing the number of subcodes N^, the 

correlation output sidelobes can be decreased significantly while providing a larger PG. 
This comes at the expense of a larger range-Doppler correlation matrix processor in the 
receiver. In addition, a longer compression time and an increase in the size of the bulk 
memory is required. 

The signal processing function used for the compression is the autocorrelation 
function (ACF). For CW radars, the compression is accomplished with a periodic 
autocorrelation function (PACF). Also, the periodic ambiguity function (PAF) is used to 
show the delay-offset versus frequency-offset in the receiver and shows the resulting 
sidelobe structure. This is especially useful for monitoring the ambiguity space for 
secondary targets that might hide in the sidelobe structure. Note that the zero Doppler cut 
of the PAF is the PACF [1]. 

C. PERIODIC AUTOCORRELATION FUNCTION 

It is important that in order to achieve high range resolution for a polyphase code 
sequence, the coded signal has to have a noise-like PACF property, close to an impulse 
function. The ACF is used for signals of finite energy. For a CW polyphase code 
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sequence with the number of phase codes , subcode duration , and a code period 
T = , a periodic complex envelope u{t) is given by [1] 

u(t) = u(t + nT) (13) 

for n = 0,±l,±2,±3,....The values of the PACF as a function of the delay r (which are 
the multiples of ) are given by 


Rirth) = ^'^uin)uin + r). (14) 

■'’c n=l 

Ideally, a perfect periodic autocorrelation function (zero-sidelobe) is desired where 
R{rti^) = l when r = 0 (mod A^), and R{rt^) = 0 when r^^O(modA^). Since the CW 
signal is continuous, a perfect periodic autocorrelation function can be achieved. 


In Figure 6a, the ACF output for a N^= 7 Barker code is plotted (in dB) with no 
noise using /^ = 1 kHz , f,= 7 kHz , cpp = 1 and A,. = 1 using the MATLAB m-file 
(amfbnV.m) [11]. The highest peak mainlobe is at the zero time delay. The width of the 
mainlobe is equal to the length of the subcode t^. For the Barker code, it is found that the 

PSL is equal to -201ogjo(A^) = -17 dB. The PACF in Figure 6b shows that the peak is 
repeated every code period T = at the value of 7. Also note from (14) that the CW 
binary Barker code is not a perfect code. 


(a) 


(b) 


Figure 6. 



Normalized correlation output of Barker code sequence N^=7 and N^=\ 
for (a) ACF and (b) PACF. 
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Today, many studies on PSK show that the PSL of the phase-eoded sequenees ean 
be redueed by ehoosing a proper phase coding. Moreover, there are CW polyphase codes 
able to achieve a zero sidelobe (perfect codes). 


D. PERIODIC AMBIGUITY FUNCTION 

In a realistic environment, the returned signal will be shifted in frequency if the 
target is moving relative to the radar’s antenna beam axis (Doppler effect). The PAF is an 
important tool to assess the response of a correlation receiver. In pulsed radar analysis, 
the ambiguity function j(r,v) is a two-dimensional function of time delay offset z and 
Doppler frequency offset v showing the correlation between the returned signal and the 
receiver matched filter. For the finite duration signal, the ambiguity function (AF) is 
determined by the properties of the pulse and the matched filter. If u{t) is the complex 
envelope of both the transmitted signal and received signal, the AF is given by [10] 


\x{t,v)\ = ^^u{t)u {t-z)e 


jlTTVt 


dt 


(15) 


where * denotes the complex conjugate. A target further from the radar than the reference 
(r = 0) position will correspond to positive z. A positive vimplies a target moving 
toward the radar. 


A more concise way of representing the ambiguity function consists of examining 
the one-dimensional zero-delay and zero-Doppler cuts; that is, ;i'(0,v) and j(r,0), 
respectively. The matched filter output as a function of a time (the signal one would 
observe in a radar system) is a delay cut, with constant frequency given by the target’s 
Doppler shift ^(r,^). Ideally, the ambiguity diagram should show the diagonal ridge 
centered at the origin and zero elsewhere, no ambiguity. 

The PAF, introduced by Levanon and Freedman [11], describes the response of a 
correlation receiver to a CW signal modulated by a periodic waveform with period 
T, when the reference signal is constructed from an integral number of periods of the 

transmitted signal (coherent processor length N^T). The target illumination time (dwell 
time) NpT must be longer than N^T. As long as the delay r is shorter than the 
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difference between the dwell time and the length of the reference signal 
0 <T <{N^^ - NJT, the illumination time can be considered infinitely long and the 

receiver response can be described by the PAF given as [8] 




1 


NX 

0 


(16) 



Figure 7. PAF of Barker code N^=l with (a) = 1 and (b) = 4. 
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In Figure 7a, the PAF of the Barker code N^=7 with N^=l is plotted using 
MATLAB m-file (amfbnV.m) from [10]. Also, the PAF of the Barker code = 7 with 
N^=4 is plotted in Figure 7b. By increasing N^, the Doppler sidelobes have been 
significantly lowered. 

In the next chapter, the PACF and PAF are used to characterize the CW Frank 
code. The residue number system is introduced and a relationship to the Frank code is 
developed. 

Since there are many variables used throughout this thesis, all of important 
variables are listed in the Appendix to provide unambiguous meanings. 
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III. FRANK POLYPHASE MODULATION AND THE RESIDUE 

NUMBER SYSTEM 

In this chapter, the Frank phase code and the residue number system (RNS) are 
introdueed. Their relationship is presented. 

A. FRANK PHASE CODE 

The Frank eode is derived from a step approximation to a linear frequeney 
modulation waveform using M frequeney steps and M samples per frequeney [4], If i 
is the sample number in a given frequency and j is the frequeney number, the phase of 
the zth sample of the 7 th frequeney is given by 

=^(<-1)0-1) (17) 

M 

where k is the time index for a single phase ehange per subeode and it ranges from 1 to 
M^,and i, j = \,2,...,M. 

The phase values and the phase values modulo 2n of the Frank eode for 
M = 8 , = 64 are shown in Figures 8 and 9, respeetively. The earrier frequeney is 

/, = 1 kHz, the sampling frequeney is /^ = 7 kHz, and the number of earrier eyeles per 
subeode is cpp =1. The ACF and the PACT are shown in Figure 10 where is the 
number of referenee eode periods in the reeeiver. The delay axis is normalized by the 
subeode period Note from Equation (14) and Figure 10b that the Frank eode is a 

perfeet eode. The PSL is PSL = 201og(l/(M;r)) or -28 dB. The PAF in Figure 11 

shows the delay-offset and Doppler-offset sidelobes. Note that the Doppler toleranee is 
refleeted in the way the sidelobes are arranged about the mainlobe. The study on Doppler 
toleranee is presented in [5] and [12]. 


19 



40 



Figure 8. Frank code modulation for M = %,N^=6A (From [1]). 



Figure 9. Signal phase (radians) modulo In versus k -index for phase change of the 

Frank code with M = 8, = 64 (From [1]). 
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(a) 



(b) 



-5 -20 

o 

o 

I -40 

o 

o 

a -60 



r/L 


Figure 10. ACF (a) and PACF (b) for the Frank code with M = 8, cpp = 1, and = 1 

(From [1]). 



Figure 11. PAF for the Fra nk code modulation with M = 8, cpp = 1, and = 1 

(From [1]). 


21 
























































B. 


THE RESIDUE NUMBER SYSTEM 


The RNS was introduced by Szabo [13]. It has been established as an important 
tool in parallel processing applications and high speed computations [14]-[16]. The RNS 
can serve as a source for extending the unambiguous range by decomposing the 
unambiguous range into a number of parallel sub-ranges (moduli) that have a smaller 
number of subcodes. Each sub-range for a different coprime modulus requires only a 
code period in accordance with that modulus. A much larger unambiguous range is 
achieved after the results of these smaller sub-range operations are recombined. Thus, by 
developing a relationship between the RNS, the polyphase distributions within the Frank 
code period, and the target detection processing, the extension of the unambiguous range 
beyond a single code period is feasible. 

The RNS residue distribution is based on a saw-tooth folding waveform, where 
the discrete residue values rise gradually and fall quickly as the input value rises 
gradually (i.e., in mod 5 the succession of discrete values for an increasing input is 0, 1, 
2, 3, 4, 0, 1, 2, 3, ...). A single RNS sequence can be generated as 

a = R^(modm) (18) 

where Q<R^<m,m is the modulus, and is the residue of a modulo m. The integer 
values within each modulus sequence are 

R„=[0,l,...,m-1,] (19) 

and are representative of the saw-tooth folding waveform. The sequence is repeated in 
both directions forming a periodic sequence with the period 

Prns = (20) 

Given a set of RNS coprime moduli, a complete residue number system is 

formed. A set of integers yQ,yj,...,y^_j(modm) form a complete residue system if 

= g(mod m) for g = 0,1,..., m -1. The dynamic range of the RNS is the number 
of paired terms from the sequence of values that contain no ambiguities and is 


22 



( 21 ) 


N 

^ RNS ~ ri 
/=1 

where m. is the /th coprime modulus and > 2 is the number of moduli, 1 < i < A^. In 
Figure 12, an RNS example for N = 2 moduli, m^=A and =5, is shown. Note the 
stair step relationship of the residues. 



input(aj 


m^= 5 



Figure 12. The RNS residue folding waveforms for m, = 4 and m^=5. 


From Equation (21), for m^=A and =5, the dynamic range is =20. In 
Table 1, the two sequences are shown along with the input a. Starting at n = 0, the total 
number of distinct paired terms is 20, 0 < a < 19. The index value Z = a + 1 in Table 1 is 
shown to count the integer number of distinct paired terms within the dynamic range. 


Table 1. Finding the RNS dynamic range for m, = 4 and = 5 . 


Input( a ) 

0 

1 

2 

3 

4 

5 

6 

7 
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10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 
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mj = 4 
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1 

2 

3 

0 

1 

2 

3 

0 

1 
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3 
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3 
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1 

2 

3 

0 

1 

2 

m2 =5 

0 

1 

2 

3 

4 

0 

1 

2 

3 

4 

0 

1 

2 

3 

4 

0 

1 

2 

3 

4 

0 

1 

2 

1 
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10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
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c. 


RESIDUE-FRANK PHASE CODE 


The RNS residues can be folded into the Frank phase code sequence by 
associating ^ with the Frank subcode phases. The phase index k starts from 1 to N^, 
where is the number of subcodes within one code period. The residue-Frank phase 
code sequence from modulo m is given by 


^m,k 


In 

m 



( 22 ) 


where k = 1,2,3,is the phase index for each phase change, and [xj indicates the 
greatest integer less than or equal to x. In Equation (22), the number of subcodes within 
one code period is N^=m^. For example, if we choose m =4, 16, the RNS folding 


waveform becomes R, 


m=4X=l,2,...,16 


= (0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3}. 


Next, the phase sequences from equations (17) and (22) are compared. For 
Equation (22), the parameters used to calculate the residue-Erank phase code sequence 
are m = 4, =16 with le = (O’ 1’ 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3}. Eorthe 

Erank code from Equation (17), the parameters are M = 4, = 16 with i, j = 1,2,3,4. In 

Table 2, the phase code sequences are shown for comparison. 

The phase code sequences of both waveforms are identical. The ACE, PACE, and 
PAE of the CW waveform generated by Equation (22) with m = 8, = 64 are also 

identical to those shown in Eigures 10 and 11, which were generated by Equation (17) 
with M = 8, A^ = 64. 
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Table 2. Comparison of the residue-Frank phase eode sequenee (m = 4) with the 

Frank phase eode sequenee {M - A, N^=\6). 


Time index 

Residue-Frank phase eode (22) 

Frank phase eode (17) 

(k) 

R 

m=4,x 

k-\ 

m 


i 

7 

.(rad) 

1 

0 

0 

^=—.0-0=0 

4 

1 

1 

27r 

^,=—0-0=0 

4 

2 

1 

0 

(!^2=—■1-0=0 
^ 4 

2 

1 

^3 =^1-0=0 

3 

2 

0 

^=—2-0=0 
^ 4 

3 

1 

«^3=^2.0=0 

4 

3 

0 

^4=^-3-0=0 

4 

1 

27r 

1^4,1 =^3-0=0 

5 

0 

1 

^=^.0-1=0 

1 

2 

2 =^0-1=0 

6 

1 

1 

«i6=--M=^ 
^4 2 

2 

2 

^2=-M=^ 

4 2 

7 

2 

1 

, 2;r 

07= - 2-1=;r 

' 4 

3 

2 

«h,2=^2-l=;r 

8 

3 

1 

^=-■3.1=^ 
^4 2 

4 

2 

^4 2 = 2^3-1=^ 

4 2 

9 

0 

2 

^=— ■0-2=0 
v-v 4 

1 

3 

^3=^0-2=0 

10 

1 

2 

o 

II 

II 

2 

3 

II 

(N 

II 

11 

2 

2 

(ill i=—-2-2=2;r 

3 

3 

«b,3=^2-2=2;r 

12 

3 

2 

^2=^-3-2=3;r 

4 

3 

^4,3=^3-2=3;r 

13 

0 

3 

'2,k 

^3=—-0-3=0 

1 

4 

27r 

(ll,4=—0-3=0 

14 

1 

3 

^4 = 2^.1-3=^ 

4 2 

2 

4 

^4 = 2^1.3=3^ 

15 

2 

3 

(ill 5 =^-2-3=3 ;r 

3 

4 

(*h,4=^2-3=3;r 

16 

3 

3 

^6=2T.3.3=9^ 

4 

4 

^4 4=2^3-3=^ 

4 2 



D. RESOLVING RANGE AMBIGUITIES USING N RESIDUE-FRANK 

PHASE CODE SEQUENCES 

I. Block Diagram of the Residue-Frank Radar System 

To extend the unambiguous detection range, N>2 coprime moduli are used with 
Equation (22) to generate the residue-Frank phase code sequences necessary to modulate 
the CW frequency. The block diagram of the radar using the residue-Frank phase code 
modulation is shown in Figure 13. In the transmitter, the first step is to generate and store 
the sequence values ^ for each modulus. Then, a direct digital synthesizer uses each 

residue sequence ^ to generate the subcodes (j)^^ ^ according to Equation (22) 

for each modulus. Each subcode corresponding to has a length of Together the 
subcodes have a code period of T. To detect the target’s unambiguous range, each 
phase code sequence must have subcodes. Each sequence is amplified and 

transmitted in succession. 


Residue-Frank 
phase code residues 


CodeWlj^r 


CodeWl2 

Code/Ki 

length 


length 

leng& 



AffliK 



•C a 

a» 

a 




Rm.k 


Direct digital 
synthesizer (22) 


Phase 


Phase 

Phase 

code W 


code 2 

code 1 


Synchronized 
codes for Bin 
detection 


Range bin 
code nil 


Range bin 
code Wj 


Rai^ebin 
code tTlif 


6 - 


o 




Amplifier 


Transmitter 


Signal 

reference 


LO 


Codel 

Compressor 


Code 2 
Compressor 


CodeN 

Compressor 


■ 6 - 

- 6 - 


ADC 


Down 

conversion 


Transmit 

antenna 


Receive 

antenna 


LNA <=> 


Figure 13. Block diagram for the radar using N residue-Frank phase code sequences. 

Upon reception, the signal is amplified and downconverted to an intermediate 
frequency compatible with an available analog-to-digital (ADC) converter technology. 
The output of the ADC is processed by a moving target indication (Elliptic) filter to 

remove the clutter, and then a phase detector is used to determine the phase and 
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magnitude of each subcode. The output of the phase detector is processed by a phase 
code compressor which is followed by a noncoherent post-detection integration (PDI) 
process. The output of the PDI is sent to a constant false alarm rate processor to detect 
and save the range bin of the target for that modulus sequence. This series of steps is 
performed for each residue-Frank phase code sequence, and the detected target’s range 
bin is sent to a RNS-to-decimal algorithm to resolve the target’s true range after all 

N residues are obtained. 

m 


2. Transmitted and Reference Codes for Compression 

In this section, the transmission of the set of N residue-Frank phase code 
sequences is discussed. The use of each phase code sequence as a reference sequence for 
code compression is discussed. An example using A = 3 moduli with m, = 3, 4, 

and m^= 5 is illustrated. The dynamic range from Equation (21) for these RNS moduli is 
Mrns = 3x4x5 = 60. 

The transmitted sequence for each modulus requires the number of subcodes 
equal to the dynamic range The transmitted signal is shown in Figure 14. For 

modulus m, =3, there are three residues within a code period (0, 1, 2). For the residue- 
Frank phase code period, the sequence has = mf = 3^ = 9. To cover a dynamic range 
of = 60 subcodes, the number of code periods required is 


^ RNS 
P,m 2 

m 


(23) 


For m, =3, N = 60/9 = 6.67 code periods. The last 3 subcodes (9x7-60 = 3) are not 
used since they are ambiguous. 

For mj = 4 and = 5, the number of subcodes for each code period are 
ml=A^ =\6 and = 5^ = 25, respectively. Also, the maximum number of transmitted 
and received code periods are 4 =3.75 (the last 4 subcodes are not used) and 
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Np 5 = 2.4 (the last 15 subcodes are not used), respectively. In Figure 14, the code period 
corresponding to each modulus is denoted as ^ where p represents code period 

number. Note that the number of received code periods processed in the range-Doppler 
matrix should be less than the number of code periods received N That is 

N ^ < N p< N Otherwise, a certain amount of correlation loss occurs. 
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Figure 14. Illustration of transmitted signal using N = 3 residue-Frank phase code 

sequences for m, = 3, mj = 4, and =5. 

In Figure 15 a, the residues for =3, =4, and m^=5 =60) are 

shown with ~^mk+ntP where n e {0,±l,±2,...}.The residue-Frank phase sequences 
from each modulus are shown in Figure 15b. Each phase code sequence has a length 
N^=m^, giving k - k+nm^- receiver uses one code period as a reference 

(N^ = 1), the reference code sequence for compression has a phase sequence from index 
= 1 to m^. 
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(a) 


T^=9 m^=3 




time-index 


(b) 

Figure 15. Plot of the residue-Frank phase eode sequenee for m, =3, m^= 4, and 

m^=5 = 60 ) showing (a) residues, and (b) phase sequenees from (22). 

In Figure 16, the eompression and range bin deteetion proeess is illustrated for 


mj = 3 and = 1. The residue-Frank phase eode sequenee is returned to the radar 

reeeiver with delay t and then proeessed by a filter matehed to a reetangular pulse of 
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duration The phase detector is used to calculate the phase of the signal which is 

passed to the compressor. The output from the compression is then processed by the 
threshold detection to see whether the amplitude reaches the desired minimum value. The 
compressor processes a single compression delay within a sampling period If the 
threshold detector detects the pulse, the presence of the target will be saved and 
associated with residue Each delay D, has a duration of (e.g. is for the time 

delay between 0 to is for the time delay between to and so on). For m^=A 

and m 3 = 5, the process is similar, except the phase reference values have to correspond 
to the transmitted phase code sequence that is being received. 


i^r — T, = Nj, = %-^\ 



* X), = fth range bin with delay (e.g. D^ is for r = 0 to is for r = <;, to 2/j,) 


Figure 16. Illustration of compression at the receiver and the range bin for RNS 

mj =3, N^= 9, and N^=l. 
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3 . 


Calculating the Target Range 


After the returned signal is eompressed, the first sample appears in the reeeiver’s 
range bin that eorresponds to a duration r. For eaeh modular sequenee, the target 
deteetion eorresponds to a residue R^. By eombining the paired residue from eaeh 
sequenee, the target’s range bin ean be ealeulated. 

Continuing the example in the previous seetion, we show the residues for 
m, = 3, m 2 = 4, and m^=5 and the eorresponding target range in Table 3. Here, the signal 

bandwidth B = \ MHz, and the subeode period t^= \ //s. The unambiguous range using 
residue-Frank phase eode is 


R 


u^RNS 


_ RNsh 
2 


(24) 


For this example, the maximum unambiguous range is R^ = 9 km = 60 ) and is 

considerably longer than that of any single code period. The range resolution is 
Ai? = / 2 = 150 m. The target’s range now has to be R< R^ in order to be detected 


with no range ambiguity. The range interval for each range bin of this example is also 
given in the Table 3. 


Table 3. Bin matrix and range intervals of 3 - channel RNS with m, =3, m^= 4, and 

m 3 = 5. 


Bin a 

m, =3 

II 

s" 

m, = 5 

Range(m) 

Bin6t 

m, =3 

3 
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II 

4^ 
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II 

s 

Range(m) 

0 
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0 
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30 

0 

2 

0 

4500^650 

1 

1 

1 

1 

150-300 

31 

1 

3 

1 

4650-4800 

2 

2 

2 

2 

300-450 

32 

2 

0 

2 

4800^950 

3 

0 

3 

3 

450-600 

33 

0 

1 

3 

4950-5100 

4 

1 

0 

4 

600-750 

34 

1 

2 

4 

5100-5250 

5 

2 

1 

0 

750-900 

35 

2 

3 

0 

5250-5400 

6 

0 

2 

1 

900-1050 

36 

0 

0 

1 

5400-5550 

7 

1 

3 

2 

1050-1200 

37 

1 

1 

2 

5550-5700 
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8 

2 

0 

3 

1200-1350 

38 

2 

2 

3 

5700-5850 

9 

0 

1 

4 

1350-1500 

39 

0 

3 

4 

5850-6000 

10 

1 

2 

0 

1500-1650 

40 

1 

0 

0 

6000-6150 

11 

2 

3 

1 

1650-1800 

41 

2 

1 

1 

6150-6300 

12 

0 

0 

2 

1800-1950 

42 

0 

2 

2 

6300-6450 

13 

1 

1 

3 

1950-2100 

43 

1 

3 

3 

6450-6600 

14 

2 

2 

4 

2100-2250 

44 

2 

0 

4 

6600-6750 

15 

0 

3 

0 

2250-2400 

45 

0 

1 

0 

6750-6900 

16 

1 

0 

1 

2400-2550 

46 

1 

2 

1 

6900-7050 

17 

2 

1 

2 

2550-2700 

47 

2 

3 

2 

7050-7200 

18 

0 

2 

3 

2700-2850 

48 

0 

0 

3 

7200-7350 

19 

1 

3 

4 

2850-3000 

49 

1 

1 

4 

7350-7500 

20 

2 

0 

0 

3000-3150 

50 

2 

2 

0 

7500-7650 

21 

0 

1 

1 

3150-3300 

51 

0 

3 

1 

7650-7800 

22 

1 

2 

2 

3300-3450 

52 

1 

0 

2 

7800-7950 

23 

2 

3 

3 

3450-3600 

53 

2 

1 

3 

7950-8100 

24 

0 

0 

4 

3600-3750 

54 

0 

2 

4 

8100-8250 

25 

1 

1 

0 

3750-3900 

55 

1 

3 

0 

8250-8400 

26 

2 

2 

1 

3900^050 

56 

2 

0 

1 

8400-8550 

27 

0 

3 

2 

4050^200 

57 

0 

1 

2 

8550-8700 

28 

1 

0 

3 

4200^350 

58 

1 

2 

3 

8700-8850 

29 

2 

1 

4 

4350^500 

59 

2 

3 

4 

8850-9000 


In Figure 17, an example of the target detection process using = 3, m 2 = 4, and 
m 3 = 5 is illustrated. The residue-Frank phase code signal is transmitted to the target and 
the target’s return is in range bin a = 5 or 770 m. In the code compression step, the 
returned signal is compressed corresponding to each modulus in consecutive order. 
Shown in each graph are residues ^ (stair-step waveform), the range bin index a , and 
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a helicopter target detection. The range bin index a goes from 0 to 59 =60). 

Since the phase code sequence is periodic, the residues is also periodic with 
Rm k = • For modulus = 3, the CC peak is repeated with = 9t^. 

The threshold detector detects the correlated signal output, and the range bin residue R^ 

is stored. Here, the range bin residue corresponding to the correlated output for channels 
mj = 3 is at R^ = 2. The processes are the same for channels 2 and 3 using the moduli 

mj = 4 and = 5, respectively. The detected target range bin residues are R^=\ and 
= 0. The residues are paired as [2 1 0], and a RNS-to-decimal algorithm is used to 
calculate the target’s true range R. From Table 3, the paired values [2 1 0] correspond to 
the range bin a = 5, and the calculated target range is between 750 m and 900 m. 


D>ii ainic 





compression time 7^ = = 25^ 


Figure 17. Example of compression output for target detection using N = 3 residue- 
Frank phase code sequences with = 3, mj = 4, and m^=5. 
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Although the look up table in Table 3 can be used to find the target’s range from 
the corresponding residues, for a practical implementation the table would be very large. 
Alternatively, the system of simultaneous congruences of the form a = R^ {modm.) can 

be solved with a straightforward application of the Chinese remainder theorem (CRT). 
The system to be solved for this example is as2(mod3), a = \{mod4), as0(mod5). 
This CRT algorithm only requires the target range bin residues . The CRT solution of 
the system is 


N l\/f 

a = Z—bA 


(25) 


where a is the range bin index. If we consider this number modulo m^, we find that 


d = —=Rm modm,.. The values can be found with the Euclidean algorithm, 
m. 

which is a repeated application of the division algorithm [17]. We note that the b. values 

are also hardwired into the system and depend only on the moduli and not at all on the 
range bin residues. 


4. Resolving Multiple Target Range Ambiguities 

In this section, the detection of multiple targets is considered. As illustrated in 
Figure 18, the residue for the two targets are periodic as ~ b^mk+nm^' there are no 

missed targets in the detection process, the sequence of residues within each modulus can 
be paired directly for the solution of the target’s range. To solve for the range of each 
target, the residues are paired as [2 2 2] for the first target and [2 0 0] for the second 
target. Using the CRT or lookup table, we find the range of the first target is 
150 < Rj < 300 m(a = 2) and the range of the second target is 3000 < Rj < 3150 m 
(a = 20). 
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Figure 18. Resolving the range ambiguity of 2 targets using RNS for 

mj = 3, m 2 = 4, and =5. 


In Figure 19, if the first set of residues for the two targets are not detected, the 
residue-Frank phase code waveform still allows the two target range values to be 
determined. Note the residues are determined but there is no information available to pair 
them up correctly. Additional information must be used in order to pair them correctly. 
For example, for m,, we have = {2, 2}, for , we have = {2, 0}, and for m^, we 

have ={2, 0}.To pair them up correctly, the magnitude information associated with 
each subcode (from the phase detector) can be used. 
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Figure 19. Missed detections of the first set of residues. 


A more convenient method is to extend or fill in the residues according to their 
code period, which can be identified through examination of the residues throughout the 
unambiguous range. This result gives the same residue configuration as shown in Figure 
18. Note that this method works even with the target residues overlapping. 
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5. Practical Considerations 

a. Detection of Target in Noise Using RNS Compression 
In this section, the performance of the N = 3 residue-Frank for m, = 3, 
= 4, and m^=5 with = 60 is examined by using MATLAB to plot the returned 
signals from two targets at ranges of 440 m and 5840 m. The transmitted signal has an 
amplitude A = 1, carrier frequency /^ = 1 MHz, sampling frequency f,=l MHz, 

cpp = 1, and ti^= \ //s. For a practical case, the signals are received with noise. The noise 

changes the returned phase sequences and affects the compression output value. Here, 
two different SNRs are examined. In Figure 20, the transmitted signal without noise for 
each channel is plotted. In Figure 21, the signal with SNR = 30 dB is plotted. The signal 
power is high compared to the noise power. Consequently, the noise signal does not 
significantly affect the phase sequences. In Figure 22, the SNR is reduced to 0 dB. The 
phase sequences are changed significantly. 
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Figure 20. Residue-Fra nk signal without noise for m, = 3, m^= 4, and m^=5 with 
^RNs = 60, A = 1, /^ = 1 MHz, f^=l MHz, cpp = 1, and t^^= \ ps. 
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Figure 21. Residue-Frank signal with SNR = 30 dB for =3, m^= 4, and m^=5 
with = 60, A = 1, /^ = 1 MFlz, /^ = 7 MFIz, cpp = 1, and t^=\ ps. 


m^= 3 



Figure 22. Residue-Frank signal with SNR = 0 dB for = 3, mj = 4, and m^= 5 with 
Mrns = 60, a = 1, /, = 1 MFIz, /^ = 7 MFlz, cpp = 1, and t^=\ ps. 
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In Figures 23 and 24, the range bin profiles for each channel after 
threshold detection of the compression output for the two targets are shown with SNR = 
30 dB and SNR = 0 dB, respectively. The first target appears at a range bin a = 2 with 
the paired residues [2 2 2] (range 300-450 m). The second target appears at range bin 
a =38 with the paired residues [2 2 3] (range 5700-5840 m). Note that the targets 
straddle the range bins since the first sample received from the targets lies in the correct 
range bin. Also shown is the noise signal that is compressed before the signals are 
returned to the receiver. Shown in Figure 24, the amplitudes of compressed pulses 
degrade when the SNR is decreased to 0 dB. Also, the sidelobes from the compression 
output are higher when the noise power is increased. Another important point is that the 
average sidelobes from the higher modulus (e.g. m^=5) are lower than those from the 

lower modulus (e.g., m, = 3 and =4). This is due to the different processing gain or 
subcodes that each sequence provides. For = 3, m^= 4, and = 5, the processing 
gains are mf =9, ml =\6, m] = 25, respectively. 
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|ACF| fromm =5 |ACF| fromm =4 |ACF| ftomm =3 



Figure 23. Target detection for SNR = 30 dB showing the range bins for m, = 3, mj = 4, and =5 with the two targets at 440 m 

and 5840 m. 
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Figure 24. Target detection for SNR = 0 dB showing the range bins for m, = 3, = 4, and =5 with two targets at 440 m and 

5840 m. 
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One situation that may occur is that the threshold detector might read the 
wrong range bin because the amplitudes of the compressed pulses fluctuate due to the 
noise present. This results in a target range error. In Figure 25, the first target in channel 
mj = 4 can be seen with the compressed pulse detected by a CFAR (threshold level at 

0.8) at residue instead of = 2. So the paired terms [2 3 2] (range 7050-7200 m) 

are computed instead of the true range at [2 2 2] (range 300^50 m). This is a large range 
error. 



Figure 25. Illustration of detection error from the threshold detector (CFAR) using 

the residue-Frank for m, = 3, = 4, and m^= 5. 


b. Required CW Signal Power Versus Required Output SNR 

In this section, the CW power required when using the residue-Frank 
phase code sequences to detect the target at the maximum unambiguous range is 
presented. Given the fixed unambiguous range = 9000 m, the CW radar power can be 

changed accordingly as a function of the SNRj^^ at the output of compressor. The 
relationship between the required SNRj^^ and the average CW power () is expressed 
as[l] 
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( 26 ) 


_ Rt{An)^k,T,F,B,,L,,L,,{SNR,J PG) 

G,G^A^ct^L^ 

where is the unambiguous range, =1.38x10^^^ joule/K (Boltzman’s constant), 
is the ambient noise temperature {T^= 290 K), is the receiver noise factor, 5^,. is the 
radar receiver’s input bandwidth in Hz, G, and are the transmitting and receiving 
antenna gains, cjj is the target’s radar cross section in m , is the two-way atmospheric 
transmission factor, is the loss between the radar’s transmitter and antenna, and 
is the loss between the radar’s antenna and receiver. 

In Figure 26, the CW signal power as a function of the required output 
signal-to-noise ratio SNR^.^ for = 3, mj = 4, and m^=5 is plotted with 5^,. = 1 MHz, 

= 1 //s, G, =G^= 30 dB, /^ = 3 GHz, /I = 0.1 m, cr^^. = 100 m^, = 5 dB, = 

= 1, and R^ = / 2 = 9000 m (maximum detection range without ambiguity). 



Figure 26. Average power of the CW transmitter for residue-Frank with 

mj = 3, mj = 4, and =5. 
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Note the difference in the average power required in each channel due to 
the different PG. The channel = 5, giving PG = N^= mf = 25, requires less power 
than the other channels, m, =3 {PG = 9) and mj =4 {PG = 16), to detect the targets 
within the same unambiguous detection range. For example, at SNRf.^ = 30 dB, the 
required CW power P^.^ for m 3 =5, m 2 =4, and mj=3 are 7, 10.5, and 18 W, 
respectively. Note that this helps contribute to the LPI characteristic of the emitter. 

c. Unambiguous Detection Range for a Constant Output SNR 

It is important to see how superior the residue-Frank phase code is 
compared to a single Frank code with the unambiguous range R^=cT / 2 = ct^N^ / 2 and 
limited to only one code period. Flere, the relationship between the average signal power 
and the maximum detection range R^^^ is 

^ ^- ‘ , 27 ) 

M / TO)L„L.. _ 

However, the resolvable detection range is limited by the unambiguous range. The radar, 
however, can detect the targets beyond R^ but there are ambiguities. 

In Figure 27, the detection range R = is plotted as a function of the 
required CW power with constant SNR^^ =13 dB, B = \ MHz, t,, = 1 jus, G^=G^= 30 dB, 
/, = 3 GHz, /I = 0.1 m, cTj. = 100 vsf, = 290 K, = 5 dB,andLj = =1- Also, 

the performance of each Frank code using the corresponding is plotted. For the Frank 
code (17), M\ = 3, N^=9 and M2 = 4, 77^ = 16,and M3 = 5, =25. Shown in the 

graph, the single Frank code has a limited unambiguous range, where Ml = 3 gives 
=3xl0'xl0 ®x9/2 = 1350 m, M2 = 4 gives = 3x10* xlO^" xl6/2 = 2400 m,and 
M3 = 5 gives R^ =3x10* xlO^® x25/2 = 3750 m. These values are less than using the 
residue-Frank phase code, giving / 2 = 9000 m. 
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Figure 27. Comparison of the maximum unambiguous range of CW radar system for 
SNRf,^ = 13 dB using the residue-Frank phase code for =3, m^= 4, and 

m 3 = 5 with the radar using each Frank code sequence individually with the 

corresponding N^. 

In conclusion, the residue-Frank phase code presented in this chapter can 
be used to extend the unambiguous range of the CW polyphase radar system. The issue 
concerned with the range bin error due to the fluctuation of the compression outputs was 
also introduced. In the next chapter, the relationship between the P4 code and the 
symmetrical number system (SNS) is presented. As will be shown in examining the PAF, 
the P4 has a better property in Doppler tolerance than the Fra nk code. This makes it more 
attractive when the emitter attempts to detect the moving targets. 
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IV. P4 POLYPHASE MODULATION AND THE SYMMETRICAL 

NUMBER SYSTEM 

In this chapter, the P4 phase eode and the SNS are introduced. Their relationship 
is presented. 

A. P4 PHASE CODE 

The P4 polyphase code is conceptually derived from a linearly frequency 
modulated waveform (LFMW). Sueh codes and compressors ean be employed to obtain 
mueh larger time-bandwidth produets. The significant advantages of P4 codes are low 
peak sidelobes and that they are more Doppler tolerant than other phase eodes derived 
from a step approximation to an LFMW [5], [12]. 

In the P4 eode the loeal oseillator frequeney, whieh is offset in the 1 and Q 
deteetors, results in a coherent double sideband deteetion. The phase sequence of a P4 
signal is given by 

(2S) 

where k is the time index for a single phase ehange per subeode and ranges from 1 to 

K- 

The phase values and phase values modulo In of a P4 eode for = 64 are 
shown in Figures 28 and 29, respeetively. The earrier frequeney is f^= 1 kHz, the 
sampling frequeney is /, = 7 kHz, and the number of earrier eyeles per subeode is cpp = 
1. The ACF and PACT for N^=l are shown in Figure 30. The peak-to-sidelobe level is 

P5L = 201ogjo ^/27(Ar^^ or -25 dB. The PACT is a perfeet code with low sidelobes. 
The PAF in Figure 31 shows the delay-offset and Doppler-offset sidelobes. 
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Figure 28. P4 phase sequence for = 64 (From [1]). 



Figure 29. Signal phase (radians) modulo In versus k -index for phase change of the 

P4 code with = 64 (From [1]). 
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Figure 30. ACF (a) and PACF (b) for the P4 eode with = 64, cpp = 1, and = 1 

(From [1]). 



Figure 31. PAF for the P4 code modulation with = 64, cpp = 1, and = l 

(From [1]). 

B. THE SYMMETRICAL NUMBER SYSTEM 

The SNS is composed of N coprime moduli. The integers within each SNS 
modulus are derived from a symmetrically folded waveform [18]. The integer values 
within each SNS modulus are derived from a mid-level quantization of the symmetrical 
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folding waveform and, therefore, ineongruent modulus m (i.e., mod 5; {0, 1, 2, 2, 1,0, 1, 
2...}). Due to the presenee of ambiguities, the set of integers within each SNS modulus 
does not form a complete residue system by themselves. It is well known that the 
inclusion of additional redundant moduli can effectively detect and correct errors that 
may occur within a RNS representation of a number. The SNS formulation is based on a 
similar concept, which allows the ambiguities that arise within the SNS to be resolved by 
using various arrangements of the SNS moduli. The vector of integer values within a 
single SNS folding waveform is given as 


and 




= 


0 ,1,..., 

m 

m 

,...,2,1,0,1 


[2 J 

L2_ 



0 ,1,...,-,--l,...,2,1,0,1 
2 2 


, when m is odd 


, when m is even 


(29) 


(30) 


where [xj indicates the greatest integer less than or equal to x. An example of the SNS 


waveforms for mj =4 and 171^=5 (A = 2)is shown in Figure 32. The period of each 
waveform equals the value of the modulus as 


^sm 


(31) 
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Figure 32. The SNS residues for m, = 4 and mj = 5. 
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The dynamic range of the SNS M is 


M = min ^ n + n ( 32 ) 

[ ^ 1=2 l=j+l j 

for coprime moduli m. , with one of the moduli even, j ranges from 1 to -1, is the 
number of moduli, and range over all permutations {2,3,..., A^}. The 

product of j from 2 to 1 is empty, and its value is 1. If all the coprime moduli in the 
system are odd, then 


M 


SNS 


= min 



1 

+ - 

2 


n 


l-j+\ 



(33) 


where j ranges from 1 to A^^-l and range over all permutations 

{1,2,3,..., A^}. Note that the dynamic range for the SNS is less than the RNS 

SNS ^ ^RNS )• 


From Equation (32), for m^=4 and = 5, the dynamic range is = 7. In 
Table 4, the two sequences are shown along with the input a. Starting at a = 0, the total 
number of distinct paired terms is 7, 0 < a < 6. The index value / = a + 1 in Table 4 is 
shown to count the integer number of the distinct paired terms within the dynamic range. 
Also shaded is the region with no ambiguity. 


Table 4. Finding the SNS dynamic range for m, = 4 and 171^=5. 


Input( a) 

0 

1 

2 

3 

4 

5 

6 

7 

8 
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11 

12 

13 
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16 
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18 

19 
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m, = 4 

0 

1 

2 

1 

0 

1 

2 

1 

0 

1 

2 

1 

0 

1 

2 

1 

0 

1 

2 

1 

0 

m2 =5 

0 

1 

2 

2 

1 

0 

1 

2 

2 

1 

0 

1 

2 

2 

1 

0 

1 

2 

2 

1 
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c. 


SYMMETRICAL RESIDUE-P4 PHASE CODE 


The SNS symmetrical residues ^ can be folded into the P4 phase sequence by 
associating 5^^ with the P4 subcode phases. The phase index k starts from 1 to N^. The 
symmetrical residue-P4 phase sequence is given by 



where k = 1,2,3, ••••, m. The number of subcodes within one code period is = m. For 
example, if we choose m = 8, A^=8, and the SNS sequence becomes 5^=8 2 g 

= {0,1,2,3,4,3,2,1}. 

Next, the phase sequences from equations (34) and (28) are compared. First, 
Equation (34) is used to calculate the symmetrical residue-P4 phase code sequence using 
m = 8, N^=S with 2 g = {0,1,2,3,4,3,2,1}. For the P4 code in Equation (28), the 

parameters are N^=8 and i = 1,2,3,....8. In Table 5, the phase sequences are shown for 
comparison. 

The phase code sequences of both waveforms are identical. The ACE, PACE, and 
PAF of the CW waveform generated by Equation (34) with m = N^=6A are also 
identical to those shown in Eigures 30 and 31 which were generated by Equation (28) 
with =64. 
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Table 5. Comparison of symmetrical residue-P4 phase code sequence 
(m = 8, N^=S) with the P4 phase code sequence {N^=S). 


Phase index 

(k) 

Symmetrical residue-P4 phase code (34) 

P4 phase code (28) 

^m=S,k 

(rad) 

k 

<t>k (rad) 
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fe = ^^2-l) _;,(2-1)=-2.7489 

0 

3 

2 

<h= 

r 

Tti 

V 2; 

2^ 

--8=-4.7124 
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g 2832 
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V 

C 8) 
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V 2 J 
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<h= 
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f ^ 
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«h=^"^^^-^(7-l)=-4.7124 
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--8=-2.7489 

4 
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2 

^^^(8-1) _;y(8-i)=-2.7489 

0 


D. RESOLVING RANGE AMBIGUITIES USING N SYMMETRICAL 

RESIDUE-P4 PHASE CODE SEQUENCES 

I. Block Diagram of the Symmetrical Residue-P4 Radar System 

The block diagram of a radar system using symmetrical residue-P4 phase code 
sequences is shown in Figure 33. This diagram is similar to the residue-Frank phase code 
radar system in Chapter III, except that now the symmetrical residue-P4 phase code 
sequences are used instead. Each phase sequence with the length equal to the dynamic 
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range is transmitted to deteet the target. In the receiver, each code sequence is 

compressed with its corresponding conjugate of phase code. 


Symmetrical Transmit 

residue-Frank phase Direct digital antenna 



Figure 33. Block diagram for the radar using N symmetrical residue-P4 phase code 

sequences. 

2. Transmitted and Reference Codes for Compression 

In this section, the transmission of N symmetrical residue-P4 phase code 
sequences is explained as well as the phase code reference for compression. An example 
using N = 3 moduli with m, = 7, mj = 8, and m^=9 is illustrated. The dynamic range, 

using Equation (32) for these SNS moduli, is = 37. 

Each sequence has the number of code length equal to the dynamic range with the 
transmitted signal as shown in Eigure 34. To create the P4 code sequence, one code 
period requires the number of subcodes = m. To cover a dynamic range of subcode, 
the number of code periods required is 

M 

• (35) 

m 

Eor m^=l and = 7, Equation (35) gives N= 31H = 5.29, which last 5 subcodes 
(7x6-37 = 5) are not included since they are ambiguous. Eor modulus m^=% and 
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mj = 9, the number of subeodes for one code period are 8 and 9 respectively. Also, 
the maximum number of transmitted and received code periods are 4.63 (the last 3 

subcodes are not used) and ^ = 4.11 (the last 8 subcodes are not used), respectively. In 
Figure 34, the code period corresponding to each modulus is denoted as ^ where p 
represents the code period number. 



Figure 34. Illustration of transmitted signal using A = 3 symmetrical residue-P4 
phase code sequences for m, = 7, m 2 = 8 , and m 3 = 9. 


In Figure 35a, the symmetrical residues 5^^ using =7, = 8 , and m 3 =9 

(Mj^j =37) are shown. Since the phase sequences are periodic, 5^^ =^mk+nm where 
n e |0,±l,±2,...}.The phase sequences from each channel are plotted in Figure 35b. Each 
symmetrical residue-P4 phase sequence has a phase code period equal to A^ = m, giving 
^mk ~‘Pm k+nm' ^hc reccivcr uses one code period (A^ = 1 ), the reference code sequence 
for compression is the phase sequence from index k = \ to m. 
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Figure 35. Plot of the symmetrieal residue-P4 for m, =l,m^= 8, and m^=9 

{Msns = 37 ) showing (a) symmetrieal residues, and (b) phase sequenees from 

(34). 


In Figure 36, the compression and range bin detection process is illustrated for 
m, = 7 and = \. The symmetrical residue-P4 phase code is returned to the radar 
receiver with delay r and then processed by a filter matched to the rectangular pulse of 
duration Then the phase detector is used to calculate the phase of the signal and pass 
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it through the compressor. The output from the compression is then processed by 
threshold detection to see whether the amplitude reaches the desired minimum value. If 
the threshold detector detects the pulse, the presence of the target will be saved with the 
symmetrical residue value S^. In and =9, the process is similar, except that 

the phase references have to correspond to the transmitted phase code that is being 
received. 



Threshold detector 
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*D^ = ith compression bin with delay 


Figure 36. Illustration of compression at the receiver and the range bin for SNS 

m, =7, N^= 7, and = \. 

3. Calculating the Target Range 

Using all channels, we see that the paired terms are mapped to a specific range bin 
input a (similar to residue-Frank code bin detection process). Continuing the example in 
the previous section for = l,m^= 8, and = 9, we show the symmetrical residue 
paired values and the corresponding ranges in Table 6. The signal bandwidth B = \ MHz, 
and the subcode period = 1 //s. The unambiguous range using symmetrical residue-P4 
is 
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( 36 ) 


R 


SNsh 


u,SNS 


For this example, the maximum unambiguous range is 5550 m =37). Range 

resolution is 150 m (AR = / 2 ). 


Table 6. Bin matrix and range intervals of 3 - channel SNS for = l,m^= 8, and 

nij = 9. 


Bin U 

s 

II 

m, =8 

3 

II 

Range (m) 

Bin U 

3 

II 

m, =8 

3 

II 

Range (m) 

0 

0 

0 

0 

0-150 

19 

2 

3 

1 

2850 - 3000 

1 

1 

1 

1 

150-300 

20 

1 

4 

2 

3000-3150 

2 

2 

2 

2 

300 - 450 

21 

0 

3 

3 

3150-3300 

3 

3 

3 

3 

450 - 600 

22 

1 

2 

4 

3300-3450 

4 

3 

4 

4 

600 - 750 

23 

2 

1 

4 

3450 - 3600 

5 

2 

3 

4 

750 - 900 

24 

3 



3600-3750 

6 

1 

2 

3 

900- 1050 

25 

3 

1 

2 

3750-3900 

7 

0 

1 

2 

1050 - 1200 

26 

2 

2 

1 

3900 - 4050 

8 

1 

0 

1 

1200- 1350 

27 

1 

3 

0 

4050-4200 

9 

2 

1 

0 

1350- 1500 

28 

0 

4 

1 

4200 - 4350 

10 

3 

2 

1 

1500- 1650 

29 

1 

3 

2 

4350-4500 

11 

3 

3 

2 

1650- 1800 

30 

2 

2 

3 

4500 - 4650 

12 

2 

4 

3 

1800- 1950 

31 

3 

1 

4 

4650-4800 

13 

1 

3 

4 

1950-2100 

32 

3 

0 

4 

4800 - 4950 

14 

0 

2 

4 

2100-2250 

33 

2 

1 

3 

4950-5100 

15 

1 

1 

3 

2250 - 2400 

34 

1 

2 

2 

5100-5250 

16 

2 

0 

2 

2400 - 2550 

35 

0 

3 

1 

5250 - 5400 

17 
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1 

1 

2550-2700 

36 

1 

4 

0 

5400 - 5550 

18 

3 

2 

0 

2700 - 2850 
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In Figure 37, an example of target detection for =7, =8, and m^=9 is 

illustrated. The target’s return is in range bin a = 4 or 650 m. Shown in the graph are 
symmetrical residue (stair-step waveform), the range bin a, and a helicopter target 

detection. The bin corresponding to the correlated output for channel m^=l is at 
=3. Note the repeated pulse from the compression output due to the periodic phase 
code. The process is the same for the channel = 8, and = 9. The residue bins are 
- 4 and = 4 respectively. Then the three symmetrical residues are paired as [3 

4 4], and the SNS-to-decimal algorithm is used to find the target’s true range R. From 
Table 6, the paired terms [3 4 4] correspond to the range bin a = 4, and the calculated 
range is between 600 m and 750 m. 



h»2 = 8 



/H3= 9 



Figure 37. Illustration of target detection by using the range bin matrix for SNS-P4 

code with m, =l,m^= 8, and = 9. 
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In Chapter III, we solved the system of eongruenees using CRT for reeovering the 
target range bin a Due to the definition of the SNS, a = ±S^ {modm^), 

a = (mod m 2 ) a = ±5^^ (modm^). For eaeh of the eongruenees, either the plus or 

minus is eorreet, but we do not know whieh. Thus, we have 2N sets of N equations. For 
example, if N = 2, we have 


i) a = 5^_(modmj), 
a = S^^(modm2), 

ii) a s5^_(modmj), 
a = -S^^(modm2), 

hi) a = -S^{modm^), 
a = -S^^(modm2), 
iv) a = -S^{modm^), 
a = 5^ (modmj). 


The CRT guarantees that eaeh of these has a unique solution modulo m^m^, and exaetly 

one of these solutions lies within the dynamie range of the system. This is the value of a. 
In faet, it is only neeessary to solve (i) and (ii) at most beeause the solutions to (iii) and 
(iv) are the negatives of the solutions to (i) and (ii), respeetively. 

Reeall that in the standard statement of the CRT we wish to solve for a where 
a = (modm.), 1 < / < A^, and the m. are pairwise relatively prime. The solution is 




^ RNS'^\'^\ I RNS'^2'^2 | | ^RNS^N^N 


m, 


m. 


m. 


(37) 


where Note that the values of h, depend only on m,. and not at all on 


/=! 


±5^ .Thus, we may assume that the eonstants c. =Mj^j^^b. / m^ are known, and the SNS- 
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to-decimal algorithm only needs to evaluate +Cj5'„,_ ±...c^5„,^ modulo and 

pieks the one value that lies within the dynamie range . 


4. Resolving Multiple Target Range Ambiguities 

In this seetion, the deteetion of multiple targets is eonsidered. In Figure 38, we 
see that the symmetrieal residue for the two targets are periodie as 5,^ ^ If there 

are no missed targets in the deteetion proeess, the sequenee of symmetrieal residues 
within eaeh modulus ean be paired direetly for the solution of the target’s range. To solve 
for the range of eaeh target, the residues are paired as [2 3 4] for the first target and [2 2 
1] for the seeond target. From the CRT or lookup table, the range of the first target is 
750 < < 900 m(a =5 ), and the range of the seeond target is 3900 < Rj ^ 4050 m 

(a = 26). 





it t 

0 12 3 3 

2 

1012332101233210123^ 

2 

1 0 1 2 3 3 2 1 0 1 


t 

it ^ 

L 

it A 

0 1 2 3 - 

3 

1101234321012343210 

1 

3 4 3 2 1 0 1 2 3 4 


t 

t a 

4 

tt A 

0 12 3- 

4 

!210123443210123443 


0 1 2 3 4 4 3 2 l|0 


Ta 

get 1 

Target 1 ; SXS [2 3 4], Bin 5 

Target 2 : SNS [2 2 1], Bin 26 

Ta 

get 2 

|0|l|2|3l4 


6 1 7 1 S 1 9 |l0|ll|l2|l3|l4|l5|l6|n|l8|l9|20|2l|22|23|24|2 


27|28|29|30|3l|32|33l34|35|36| 



\7 



Rx 


-=8 


mi=9 


Figure 38. Resolving the range ambiguity of 2 targets using SNS for 

mj =1 ,m^= 8 , and m 3 = 9. 


In Figure 39, if the first set of symmetrieal residues for the two targets are not 
deteeted, the symmetrieal residue-P4 phase eode still allow the two target range values to 
be determined. Note the symmetrieal residues are determined but there is no information 
available to pair them up eorreetly. Additional information must be used. For example, 
for m,, we have Sj={2,2}, for m^, we have 5g={3,2}, and for m^, we have 
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= {4, 1}. To pair them up correctly, the magnitude information associated with each 
subcode (from the phase detector) can be used. 

Rx 

m 1=7 

m2=8 


Bin a 

Figure 39. Missed detections of the first set of symmetrical residues. 

A more convenient method is to extend or fill in the symmetrical residues 
according to their code periods, which can be identified through examination of the 
symmetrical residues throughout the unambiguous range. This result gives the same 
symmetrical residue configuration as shown in Figure 38. Note that this method works 
even with the target symmetrical residues overlapping. 

5. Practical Considerations 

a. Detection of Target in Noise Using SNS Compression 

In this section, the performances of the N = 3 SNS for m, =1, m 2 = 8 , and 

m^=9 with MSNS ~ examined by using MATLAB to plot the returned signal from 

two targets at ranges of 890 m and 2990 m, respectively. The transmitted signal have 
amplitude A = 1,carrier frequency /, = 1 MFlz,sampling frequency f,=l MHz, cpp = 1, 

and = 1 //s. The two different SNRs are examined. In Figure 40, the transmitted signal 
without noise for each channel is plotted. In Figure 41, the signal with SNR = 30 dB is 
plotted. For the data shown in Figure 42, the SNR is reduced to 0 dB. The phase 
sequences are changed significantly with the increase in noise power. 
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Figure 40. Symmetrical residue-P4 signal without noise for m, = l,m^= 8 , and 
m 3 = 9 with = 37, A = 1, /^ = 1 MHz, /^ = 7 MHz, cpp = 1, and = 1 //s. 
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Figure 41. Symmetrical residue-P4 signal with SNR = 30 dB for -1 8 , and 

m 3 = 9 with M^ns = 37, a = 1, /^ = 1 MHz, f,=l MHz, cpp = 1, and t^- \ //s. 
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Figure 42. Symmetrical residue-P4 signal with SNR = 0 dB for =1, m 2 = 8 , and 
m 3 = 9 with = 37, A = 1, /^ = 1 MFlz, /^ = 7 MHz, cpp = 1, and = 1 ps. 


In Figures 43 and 44, the range bin profiles for each channel are shown 
using SNR = 30 dB and SNR = 0 dB, respectively. The first target appears at range bin 
a = 5 with the paired symmetrical residues [2 3 4] (range 750-900 m). The second target 
appears at range bin a = 19 with the paired symmetrical residues [2 3 1] (range 2850- 
3000 m). Also shown is the noise signal that is compressed before the returned signals 
from the targets arrive at the receiver. In Figure 44, the amplitudes of compressed pulses 
degrade when the SNR is decreased to 0 dB. The effect from the noise signal degrades 
the ability to detect the targets similar to that discussed in the residue-Frank phase code 
radars. The average sidelobes from the higher modulus (e.g. m^=9) are lower than those 

from the lower modulus (e.g., m^=7 and mj = 8 ) due to the different processing gain. 
For m, = 7, mj = 8 , and = 9, the processing gains are = 7, = 8 , and = 9, 

respectively. 
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|ACF| fnomm =9 |ACF| fmmm =8 |ACF| fromm = 


3 3 



Figure 44. Target detection for SNR = 0 dB showing the range bins for = l,m^= 8,and m^=9 with two targets at 890 m and 

2990 m. 
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Like the RNS phase eode, a situation in where the threshold deteetor 
might read the wrong residue bin ean oeeur beeause the amplitudes of the eompressed 
pulses fluetuate due to the noise. This results in the target range error. In Figure 45, the 
first target in ehannel mj = 8 is seen with the eompressed pulse deteeted by CFAR 

(threshold level at 0.85) at symmetrieal residue S^=2 instead of =3. So, the paired 

terms are now [2 2 4] (no range mateh) instead of the true range at [2 3 4] (range 750- 
900 m). This results in target range error. 



Figure 45. Illustration of deteetion error from the threshold deteetor (CFAR) using 
the symmetrieal residue-P4 for m, = l,m^= 8, and = 9. 

b. Required CW Signal Power Versus Required Output SNR 

In this seetion, the CW power required when using the symmetrieal 
residue-P4 phase eode sequences to detect the target at the fixed unambiguous range 
R,, = 5550 m is presented. The relationship between the required SNRj^^ and the average 

CW power is given by Equation (26). 

In Figure 46, the CW signal power as a function of the required output 
SNRf,^ for mj =7, m^= 8, and m^=9 is plotted with 5^,. = 1 MHz, = 1 jUS, 
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G,=G^= 30 dB, = 3 GHz, A = 0.lm,a^= 100 m^ F^=5 dB, = 1, and 

11 = 5550 m (maximum detection range without ambiguity). 



Figure 46. Average power of the CW transmitter for symmetrical residue-P4 with 

nij =1 ,m^= 8 , and m 3 = 9. 

Note that the difference of the average power used in each channel is due 
to the PG. The channel = 9, giving PG = N^= m. = 9, requires less power than the 
other channels, m^=l {PG = 7) and mj = 8 (PG = 8). 

c. Unambiguous Detection Range for a Constant Output SNR 

For a single P4 code, the unambiguous range is limited to only one code 
period = cT / 2 = ct^N^ /2). In this section, we compare the maximum unambiguous 
range using the symmetrical residue-P4 phase code sequences with the corresponding 
single P4 code sequences. The relationship between the average signal power P^^ and 

the maximum detection range is given by Equation (27). 
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In Figure 47, the detection range is plotted as a function of the required 
CW power using N = 3 SNS for m^=7, m^=%, and m.^=9 with SNR^^=\3 6E, 

B = 1 MHz, t, = 1 //s, G,=G^= 30 dB, /^ = 3 GHz, A = 0.lm,cx^ =100 m\T^= 290 K, 
F^=5 dB,and - ^rr comparison, the performance of each P4 code 

using the corresponding is plotted. For the P4 code from Equation (28), = 1 and 

~ 8 , and = 9. The single P4 code has a limited unambiguous range, = 7 gives 

=3x10^x10 ®x7/2 = 1050 m, 77 ^ 2-8 gives i?,, = 3x10*^ xlO *^ x8/2 = 1200 m, and 

77 ^ 3=9 gives =3x10* x 10x9/2 = 1350 m. These values are less than using the 

77 = 3 symmetrical residue-P4 phase code sequences, giving / 2 

= 5550 m. 



Figure 47. Comparison of the maximum unambiguous range of a CW radar system 
for SNRjf^ = 13 dB using symmetrical residue-P4 phase code for = 7, 

= 8, and m^=9 with the radar using each P4 code individually with the 

corresponding 77^. 
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In conclusion, the symmetrical-P4 phase code sequence presented in this chapter 
can be used to extend the unambiguous range of the CW polyphase radar systems. The 
issue concerned with the range bin error due to the fluctuation of the compression outputs 
was considered. In the next chapter, the relationship between the P4 code and the robust 
symmetrical number system (RSNS) is presented to examine the integer Gray code 
properties and how they can prevent an error in the detected symmetrical residues. 
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V. P4 POLYPHASE MODULATION AND THE ROBUST 
SYMMETRICAL NUMBER SYSTEM 

As discussed in Chapters III and IV, using the RNS and the SNS may give a large 
error in the target’s range when the threshold detector detects the wrong residue due to 
the presence of noise which can change the phase of the received signal. In this chapter, 
the RSNS is presented, and it is shown that the range errors in the robust symmetrical 
residue-P4 phase code sequence can be bounded to a small value due to the Gray code 
property of the RSNS [20]. The P4 phase expression is presented and examples are given 
to illustrate the target detection process with this new waveform. 

A. P4 PHASE CODE 

The details of the P4 phase code are discussed in Chapter IV. The P4 phase code 
is derived from a linearly frequency modulated waveform (LFMW). To calculate the 
phase sequence. Equation (28) is used. Recall the characteristics of the P4 code sequence 
is its superior PAF sidelobe structure and its Doppler tolerance properties. 

B. THE ROBUST SYMMETRICAL NUMBER SYSTEM 

The RSNS is a modular system consisting of A > 2 integer sequences with each 
sequence associated with a coprime modulus m.. Due to the presence of ambiguities, the 
set of integers within each RSNS sequence do not form a complete residue system. A set 
of integers yo,yj,...,y^_j(modm) form a complete residue system if = g(modm) for 

g = 0,l,...,m-l [21]. The ambiguities within each modulus sequence are resolved by the 
use of additional moduli and considering the vector of paired integers from all N 
sequences. The RSNS is based on the sequence 

[RSlt] = [0,\X---,tn-l,m,m-\,...,2,\'\ (38) 

for an integer h such that 0 < h < 2m. To form the A-sequenee RSNS, each term in (38) 
is repeated A times in succession. The integers within one folding period of a sequence 
are then 
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( 39 ) 


m, m, • • • m, m -1, • • •, m -1, • • •, 1, • • •, 1], 

This results in a periodic sequence with a period of [20] 

= 2Nm. (40) 

That is, each RSNS sequence contains m. integers with the folding period of each 
sequence equal to 2Nm-. The construction of the N sequences ensures that any two 

successive RSNS vectors (paired integers from the N sequences) differ by only one 
integer, resulting in an integer Gray code property. 

Each sequence corresponding to m. is right (or left) shifted by s. =i-\ places. 
The chosen shift values {, • • •, ■s'a? ) must also form a complete residue system modulo 
N. Further, if each sequence is extended periodically with period 2Nm as 
m,h whcrc u G {0, ±1, ±2,...}, thcn RS^J^ is called a symmetrical residue of 
{h + n2Nm ) modulo 2Nm. 

The calculation of and its location in the sequence are a function of N and 

the choices for m. and s-. A closed-form solution for computing for N = 2 

moduli is reported in [22], and summarized as follows: 


For m, > 3 and mj = m, +1 



^RSNS ~ 3(mj +m2)-6 = 6m^ -3. 

(41) 

For 5 < m, 

< and m 2 < m, + 2 



^RSNS + 2m2 — 5. 

(42) 

For 5 < m, 

< and m 2 > m, + 3 



^RSNS ^ 4nij + 2m2 — 2. 

(43) 


A closed-form solution for A = 3 moduli of the form 2" -1,2", 2"^ +1 is 

3 2 15 „ .... 

^RSNS^^^l (44) 
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where z is any integer and m, > 3 [23]. An effieient algorithm for eomputing and 

its position for any general set of moduli is reported in [24]-[25]. 

The shift values y for eaeh sequenee do not affeet the RSNS dynamie range but 
do make a differenee in the loeation of the beginning position h of the dynamie range 
and the ending position h + -1. After ehoosing the beginning h , it is useful to 

know the symmetrieal residues at this point (e.g., to align the symmetrieal residues for 
eaeh sequenee). Sinee eaeh folding period eonsists of 2Nm^ integers, the symmetrieal 
residues are determined by first subtraeting off an integer number of 2Nm^ integers as 


n, =h — 


2Nm. 


2Nm. . 


(45) 


This value is then used to find the symmetrieal residue RS^ ^ as [20] 




n, - S; 


N 


2Am; + N -n. + Si -I 
N 


Si < til - + ■^, + 1 

Am, + S:+ 2 < n. < 2Am, + 5',. -1. 


(46) 


Let be the veetor of A paired integers from eaeh sequenee at h where h is the 
position of the veetor X,^ within the RSNS. In Table 7, h and X^^ for an A = 3 RSNS 
system with m^ =[3 4 5]^ and right shift Si =[0,1,2]^ are shown. For this example, the 
RSNS dynamie range =43, and the position begins at h = 6l [20]. The veetor 

Agj =[2 4 1]^. The set of integers that lie within the dynamie range =43 eontain 

no ambiguities. Also note in Table 7 the integer Gray eode property where any eode 
transition results in just one integer ehanging value by ±1. Also shown is the bin index a 
whieh runs from 0 to 42 (total 43) that will be used to refer to the range bin deteetion. 
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Table 7. The RSNS folding waveforms for m^=3 ( 5 j = 0), m^=A {s^ = 1), and 

m^=5 (53 = 2) (After [20]). 
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C. ROBUST SYMMETRICAL RESIDUE-P4 PHASE CODE 

Similar to the SNS phase eode, the RSNS folding waveform ean be related to the 
P4 eode sequenee. The RSNS symmetrieal residues RS^^. ean be folded into the P4 

phase sequenee by assoeiating RS^ with the P4 subeode phases. The phase index k' 
starts from 1 to A , where N is the number of RSNS subeodes within one eode 

^RSNS ^RSNS 

period. The robust symmetrieal residue-P4 phase sequenee is given by 


^m.k' ~ 


n 

2 m 




n 

-m 

2 


(47) 


and 


A, = 2mA 


(48) 


where A is the number of ehannels. For example, if we ehoose m = 4, A = 2, and 
=16. The robust symmetrieal residue-P4 phase eode sequenee is ^5^=4*..= 
{0,0,1,1,2,2,3,3,4,4,3,3,2,2,1,1}. Note that the folding sequenee is similar to SNS with 
= = { 0 , 1 ,2,3,4,3, 2 , 1 } exeept that the robust symmetrical residue-P4 phase is 

repeated A times; there are A phase values for each robust symmetrical-P4 phase for 
every P4 phase value. That is, both k' = \ and k' = 2 of the RSNS corresponds to A: = 1 
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for the P4 phase. Consequently, the number of subeodes for the P4 eode is related to 
the number of RSNS subeodes as 

^RSNS 

N (49) 

" N 

Also, the PG for the robust symmetrieal residue-P4 waveform is 

= (50) 

1. Example of Phase Sequence Calculation 

In this section, the phase code sequences from equations (47) and (28) are 
compared. Equation (47) is used to calculate the phase sequence for m = 4, N = 2, and 

^‘:rsns Using the RSNS, we find that the robust symmetrical residues are 

PSm= 4 k' “ {0^ 0, 1, 1, 2, 2, 3, 3, 4, 4, 3, 3, 2, 2, 1, 1} where k' = \,2,...,\6. In order to 

compare to the P4 phase sequence in Equation (28) with IN = \6/2 = S, k = 

1,2,3,...,8, the two phase sequences are shown in Table 8. 


Table 8. Comparison of the robust symmetrical residue-P4 phase code sequence 

(m = 4, A = 2, A =16) with the P4 phase code (N= 8). 

\ ^RSNS ^ \ C 


Phase 

index 

(U) 

Robust symmetrical reside-P4 phase code 

(47) 

P4 phase code (28) 

P^m=4,k' 

^m=4,r(rad) 

k 

A (rad) 

1 

0 


1 

o 

2 

0 


3 

1 


2 

2 

^ = ^U-l) ^(2 1)= 2.7489 

8 

4 

1 
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5 

2 

<h=[ 


--4=-4.7124 
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2. ACF, PACF, and PAF 

In this section, we examine the ACF, PACF, and PAF of the signal. To compare 
the robust symmetrical reside-P4 phase code sequence with the P4 phase code sequence, 
the ACF, PACF, and PAF are examined for each waveform. In Figure 48, the robust 
symmetrical residues using m. =16, A = 2, and =64 are shown. In Figure 49, the 

robust symmetrical residue-P4 phase code sequence residues are shown (using Equation 
(47)) with /^ = 1 kHz, f^=l kHz, and cpp = 2. Note that cpp represents the number of 
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carrier cycles per phase change, and the cpp value for the P4 is N times that of the 
robust symmetrical residue-P4 phase code. For the P4, the subcode period is 



For the robust symmetrical residue-P4, the subcode period is 


(51) 



(52) 


From (51) and (52), we have 


h=Ntl. (53) 

Using Equation (48), we see that the phase sequence from the robust symmetrical 
residue-P4 for m = \6, N = 2,and N = 64 is actually the same as the P4 phase code 

' ‘^RSNS ■' ^ 


sequence using = 32. 



Figure 48. RSNS residues RS^ j. for modulus m = 16, N = 2, and = 64. 

C m,K y y (^RSNS 
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Figure 49. The robust symmetrical residue-P4 phase code sequence using m = 16, 

N = 2, ~ shown with the relationship with the P4 code using = 32. 

The power spectral density (PSD) of the signal is shown in Figure 50. The signal 
bandwidth B is reduced by the factor of N. With cpp = 2, the bandwidth of the robust 
symmetrical-P4 signal is 5 = /,/cpp = 1/6, =500 Hz (instead of 5 = 1 kHz when 
cpp =1). In Figure 51, the ACF and PACF are plotted. The signal’s characteristics are 
the same as the P4 phase code sequence using = 32, t^=2 ms. In the ACF plot, the 

output gives PSL = 201ogjo |.^2/ (A^;r^) j = -22 dB. Note that the width of the mainlobe 

in the ACF plot is 2t^ and is proportional to the number of moduli N. That is, the higher 
the number of moduli A, the broader the mainlobe. In Figure 52, the PAF is plotted. The 
PAF from the robust symmetrical residues using m.=\6,N = 2, and =64 in 

Figure 52a is compared to the PAF from the P4 code sequence using =32 in Figure 
52b. Both PAFs are the same. 
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Figure 50. 
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Figure 51. The ACF and PACF of RSNS phase-coded signal 
(m = \6, N = 2, N = 64 ) at the receiver with A = 1, 

^ ^RSNS ^ ^ 

/^ = 1 kHz, /j = 7 kHz, and cpp = 2. 
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(a) 



(b) 

Figure 52. The PAF of the signal for (a) the robust symmetrieal residue-P4 phase 
code sequence using m = 16, N = 2, = 64, and tl=l ms, and (b) the P4 code 

using = 32, and t^=2 ms. 
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D. RESOLVING RANGE AMBIGUITIES USING N ROBUST 

SYMMETRICAL RESIDUE-P4 PHASE CODE SEQUENCES 

I. Block Diagram of the Robust Symmetrical Residue-P4 Radar System 

To extend the unambiguous detection range, N>2 coprime RSNS sequences are 
used. The block diagram of the radar using V - channel RSNS phase coding is shown in 
Figure 53. This diagram is similar to the residue-Frank phase code radar system except 
that now the RSNS moduli are used instead. Each phase sequence with the length equal 
to the dynamic range M is transmitted to detect the target. In the receiver, each code 
sequence is compressed with its corresponding reference phase code. 
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Figure 53. Block diagram for the radar using N robust symmetrical residue-P4 phase 

code sequences. 


2. Transmitted and Reference Codes for Compression 

In this section, the transmission of a set of N robust symmetrical residue-P4 
phase coded sequences is explained. Also, the use of the phase code sequences as the 
reference sequence for code compression is discussed. An example using A = 3 RSNS 
moduli with m, = 3, = 4, and m^=5 is illustrated. The dynamic range for these RSNS 

moduli from Equation (44) is =43. The starting position h and ending position 


79 

























































h + -1 are 61 and 103, respectively [20]. The robust symmetrical residues within 

the RSNS moduli are paired as shown in Table 7. 

The transmitted signal is shown in Figure 54. Each sequence has a length equal to 
the dynamic range. For = 3, there are four robust symmetrical residues within the 

sequence, RS^ =0, 1, 2, and 3. To create the P4 code sequence, one code period needs 
the number of subscodes N =2Nm, =2x3x3 = 18. For m, =4, the robust 

symmetrical residues are 0, 1, 2, 3, and 4. The required number of subcodes is 
2x3x4 = 24. For =5, the robust symmetrical residues are 0, 1, 2, 3, 4, and 5. The 
number of subcodes is 2x3x5 =30. As shown in the Figure 54, the beginning of each 
sequence is not necessarily the beginning of the phase code sequence k'= \ (e.g., the 
beginning phase sequences for m, =3, m 2 = 4 , and m^=5 are k' = S,k' = \3, and 
k' = 30, respectively. The reference code sequences in the compressors can be chosen to 
be either starting at the beginning of the sending signal waveforms or the starting of 
phase sequence at k' = 1. The importance is that the reference code has to have a length 
equal to the code period for each channel in order to correctly calculate the code 
compression output. In this example, we choose to use the reference signals starting from 
the beginning of the transmitted signal (^4 to (j)^ for = 3, to for = 4, and 

(/)^Q to ^2^29 for mj = 5). The number of sequence periods is calculated as 

^ ^RSNS _ ^RSNS 

N 2mN ■ ^ ^ 

^RSNS 

Sending 43 subcodes, the modulus m^=3 needs 3 = 43 /18 = 2.4. For modulus 
m 2 = 4 and m 3 = 5, the number of subcodes are 24 and 30, respectively. Also, the 
number of transmitted code periods are =1.8 and =1.4, respectively. Here, 
only channel m, = 3 is sending more than one complete code period. So we can expect to 
see two compressed pulses from channel m, = 3 if all subcodes return back to the 
receiver (there will be only one pulse in the other two channels). 


80 



= 3 


Notused 


Wo = 4 


Wj = 5 


RS,j,. = 


RS,^. ■■ 
k'- 

A.= 


RS, 


0 1 2 3 4 5 

6 7 8 9 10 11 12 13 14 15 16 1" 18 

19 

20 

21 

:: 

23 

24 

2} 

26 


28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 ■ 

|2 2 3 3 3 2 

2 2 1 1 l|0 0 0 1 1 1 2 2 

2 

3 

3 

3 

2 

2 

2 

1 

1 

i| 

» 

0 

0 

1 

1 

1 

2 

2 

2 

3 

3 

3 


8 9 10 11 12 13 

14 15 16 17 18 1 2 3 4 5 6 ' 8 

9 

10 

11 

12 

13 

14 

15 

16 

r 

18 

1 

2 

3 

4 

5 

6 


8 

9 

10 

11 

12 

13 14 

y4k ^ A ^\y 

















1 








Tj=m\ 








1 

A 

2 















0 1 2 3 4 5 

6 7 8 9 10 11 12 13 14 15 16 17 18 

19 

20 

21 


23 

24 


26 

2- 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


|4 4 4 3 3 3 

2 2 2 1 1 l| 0 0 0 1 I I 2 

2 

2 

3 

3 

3 

4 

4 

4 

3 

3 

3 

2 

2 

2 

1 

I 

i| 

[o 

0 

0 

1 

I 

>■ 

13 14 15 16 r IS 

19 20 21 22 23 24 1 2 3 4 5 6 * 

8 

9 

10 

11 

12 

13 

14 

j5 

16 

1? 

18 

19 

20 

21 


23 

24 

i 

2 

3 

4 

5 

6 * 

^ 4h ^ ^ 4{[s4{\4{'^ 




















^.1 “ * 
























0 1 2 3 4 5 

6 ' 8 9 10 11 12 13 14 15 16 17 18 

19 

20 

21 


23 

24 

25 

26 

2- 

28 

29 

30 

31 

32 

33 

34 

35 

36 

3" 

38 

39 

40 

41 ■ 

|l|0 0011 

1 2 2 2 3 3 3 4 4 4 5 5 5 

4 

4 

4 

3 

3 

3 

2 

2 

2 

I 

I 


a 

0 

0 

I 

I 

1 

2 

2 

2 

3 



Notused 


30 i : 3 4 5 6 ' t 9 10 u i: 13 :4 15 :< P n 19 xt :i :: 23 24 

^^4o^ ^4 A ^ ^ ^4 4 

T,., = 30l\ 


2t 29 30 1 
!? 


“ '• Notused 


Figure 54. 


Illustration of transmitted signal using N = 3 robust symmetrieal residue- 
P4 phase eode sequenees for m, = 3, mj = 4, and m^= 5. 


In Figure 55a, the RSNS folding waveforms for m, = 3, m^= 4, and m^=5 with 
^RSNs ~ plotted. Also, the phase sequenees from eaeh channel are plotted in 

Figure 55b. Since RS,„j^^„ 2 Nm = ^^m,h where n g{0 ,±l,±2,...},the phase sequences are 
periodic as = <Pn,,k’ • 


In Figure 56, the compression and range bin detection process is illustrated for 
RSNS mj =3, and N^=l. The returned phase coded signal is first processed by a filter 

matched to the rectangular pulse of duration The phase detector is used to calculate 
the phase of each sample and then pass it through the compressor. The output from the 
compression is then processed by the threshold detection to see whether the amplitude 
reaches the desired minimum value. If the threshold detector detects the pulse, the 
presence of the target will be saved with the robust symmetrical residue value R5,„.In 

1112=4 and =5, the process is similar, except that the phase references have to match 
with the transmitted phase code that is being received. 
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Figure 55. Plot of the RSNS using =3, m^= 4, and m^=5 with M~ bor 
(a) folding waveforms, and (b) their phase sequences. 
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Figure 56. Illustration of compression at the receiver and the range bin for RSNS 

m, =3, = 18, and N^=l. 


3. Calculating the Target Range 

Using all sequences, we map the paired terms to a specific range bin a (similar to 
the residue-Frank code bin process in Chapter III). Continuing the example in the 
previous section for m, = 3, mj = 4, and = 5, we show the robust symmetrical residue 
paired values and the corresponding ranges in Table 9. The signal bandwidth is 
5 = 1 MHz, and the subcode period is N = 0.33 //s. The new expression formula 

for calculating the unambiguous range using RSNS is 


R. 


cM t' 

RSNS'-b 


II, RSNS 


(55) 


For this example, the maximum unambiguous range is 5^ ^^^^ = 2150 m. The range 
resolution is 50 m (AR = / 2 ). 
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Table 9. Bin matrix and range intervals of 3 - channel RSNS for = 3, mj = 4, and 

mj = 5. 


Bin U 

m, = 3 

nij - A 

m^=5 

Range (m) 

Bin d 

m, =3 

niy =4 

m^=5 

Range (m) 
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In Figure 57, an example of target detection for = 3, mj = 4, and m^=5 is 
illustrated. The return signal is in the range bin a = 5 or 270 m. First, the phase code 
signal of the modulus m, = 3 (top) is compressed. Recall that the RSNS has a width of 
the compression output that covers N = 3 code sequences at robust symmetrical residues 
RSj^ = 3, 2, and 2. So the threshold detection has to decide to take one bin, e.g., at the 
peak of the compressed pulse, giving the robust symmetrical residue = 2. The 
process is the same for the channels mj = 4 and =5. The robust symmetrical residues 
are ^^4 = 3 and RS^ = 1, respectively. The three robust symmetrical residues are paired 

as [2 3 1] and must be converted to the decimal RSNS bin number. In Table 9, the 
corresponding range is between 250 m and 300 m (bin a = 5). 





Figure 57. Illustration of target detection by using the range bin matrix for RSNS-P4 

code with m, = 3, mj = 4, and =5. 
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4. 


Resolving Multiple Target Range Ambiguities 


In Figure 58 for =3, the compressed pulses are repeated as the returned code 

runs through the reference code in the compressor. The two target returns overlap one 
another at range bins a = 19, 20, and 21 and are ambiguous. Next we resolve the range 
ambiguity by using 3 channels. If the threshold detector is set to measure at the peak of 
the compression output, the robust symmetrical residues from the first target are paired as 
[3 4 0], giving a range bin a = 2 (range 50-100 m). Similarly, the robust symmetrical 
residue paired terms from the second target are [3 2 4], giving a range bin a = 20 (range 
1,000-1,050 m). The two targets have different paired terms, so the range ambiguity can 
be resolved. 
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Figure 58. Resolving the range ambiguity of 2 targets using RSNS for 

mj = 3, mj = 4, and =5. 


As mentioned in Chapters III and IV, we can resolve the targets even if the first 
set of the residues are not detected. Imagine we use a higher set of RSNS moduli and the 
dynamic range is higher with the result that the code period repeats. To find the true 
range of the targets that are missed, the magnitude information or the filling in the range 
bin according to the code periodicity can be used. 
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5. Practical Considerations 

a. Detection of Target in Noise Using RSNS Compression 

In this section, the performances of the N = 3 RSNS for m, = 3, mj = 4, 

and m^=5 with Mj^^ns ~ examined by using MATLAB to plot the returned signal 

from two targets at ranges of 140 m and 1040 m. The transmitted signal has amplitude 
A = 1,carrier frequency MHz, sampling frequency f,=9 MHz, cpp = 1, = 1 ps, 

and tl =0.33 ps. In Figure 59, the transmitted signal without noise for each channel is 

plotted. In Figure 60, the signal with SNR = 30 dB is plotted. The signal power is high 
compare to the noise power. As a result, the noise signal does not significantly affect the 
phase sequences. In Figure 61, the SNR is reduced to 0 dB. The phase sequences are 
changed significantly. 
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Figure 59. 
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Robust symmetrical residue-P4 signal without noise for 
= 4, and m^=5 with =43, A = 1, f =1 MHz, /, = 9 MHz, 
cpp = 1, and tl = 0.33 ps. 
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Figure 60. Robust symmetrical residue-P4 signal with SNR = 30 dB for 
m, = 3, mj = 4, and m^=5 with = 43, A = 1, /. = 1 MHz, /, = 9 MHz, 

cpp = 1, and tl = 0.33 //s. 
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Figure 61. Robust symmetrical residue-P4 signal with SNR = 0 dB for 
m, =3, m 2 = 4, and mj = 5 with = 43, A = 1, = 1 MHz, f, =9 MHz, 

cpp = 1, and tl = 0.33 ps. 
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In Figures 62 and 63, the range bin profiles from two targets are shown 
using SNR = 30 dB and SNR = 0 dB, respeetively. The first target appears at range bin 
a = 2 with the paired robust symmetrieal residues [3 4 0] (range 100-150 m). The seeond 
target appears at range bin a = 20, giving paired terms [3 2 4] (range 1000-1050 m). As 
shown in Figure 63, the amplitudes of the eompressed pulses degrade when the SNR is 
deereased to 0 dB. Also, the sidelobes from the eompression output are higher when the 
noise power is inereased. Another important point is that the average sidelobes from the 
higher modulus (e.g. m^=5) are lower than those from the lower modulus (e.g. m^=3 

and m 2 =4) due to the different proeessing gains. For =3, =4, and =5, the 

proeessing gains are = 2mj = 6, A ^2 = = 8, and N^^=2m^=\0, respeetively. 


89 



|ACF| from m =5 |ACF| from m =4 |ACF| ftom m =3 
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Figure 62. Target detection for SNR = 30 dB showing the range bins for = 3, m^= 4, and m^=5 with two targets at 140 m and 

1040 m. 
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|ACF| from m = 5 |ACF| from m = 4 |ACF| from m. 
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Figure 63. Target detection for SNR = 0 dB showing the range bins for = 3, m^= 4, and m^=5 with two targets at 140 m and 

1040 m. 
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Using N robust symmetrical residue-P4 phase code sequences allows 
control of the range coding error. The threshold detector might read the wrong robust 
symmetrical residue bin because the amplitudes of the compressed pulses are fluctuating. 
However, the paired terms from the robust symmetrical residues will change only one 
position (the Gray code property). In Figure 64, the compressed pulse from the first target 
in mj = 4 is detected by CFAR (threshold level at 0.85) at the robust symmetrical residue 

^^4 = 3 instead of ^^4 = 4. The paired terms are now [3 3 0] (range 150-200 m) instead 

of the true range at [3 4 0] (range 100-150 m). The maximum range error is only one 
range bin or 50 m. 



Figure 64. Illustration of detection error from the threshold detector (CFAR) using 
the robust symmetrical residue-P4 for m^ =3, m^= 4, and =5. 


b. Plot of Required Signal Power Versus Required SNR Output 

In this section, the power required when using the robust symmetrical 
residue-P4 phase code sequence to detect the target at the maximum unambiguous range 
is presented (given the fixed unambiguous range R„=2150m). The relationship 

between the required SNRj^^ and the average CW power P^^ is given by Equation (26). 

In Figure 65, the CW signal power as a function of the required SNRj^^ for 
m, = 3, m 2 = 4, and m^=5 is plotted with 5^,. = 1 MHz, C = 1 jUS, t[ = 0.333 //s. 
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G,=G^= 30 dB, /^ = 3 GHz, A = 0.lm,a^= 100 m^ F^=5 dB,and = 1, 

and = Ru - rsns^I / 2 = 2150 m (maximum detection range without ambiguity). 



Figure 65. Average power of the CW transmitter for robust symmetrical residue-P4 

with mj = 3, mj = 4, and m^=5. 

Note that the difference of the average power used in each channel is due 
to the PG. The channel = 5, giving PG = 2m = 10, requires less power than the other 

channels, m, = 3 ( PG = 6) and mj = 4 ( PG = 8 ). 

c. Unambiguous Detection Range for a Constant Output SNR 

For a single P4 code, the unambiguous range is limited to only one code 
period duration = cT / 2 = / 2. In this section, we compare the maximum 

unambiguous range using the robust symmetrical residue-P4 phase code sequence with 
the corresponding single P4 code sequences. The relationship between the average signal 
power P^^ and the maximum detection range R^^^^ is given by Equation (27). 
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In Figure 66, the detection range is plotted as a function of the required 
CW power using N = 3 RSNS for = 3, m^= 4, and m^=5 with a constant 

=\3dB, 5 = 1 MHz, t, =l//s, t; = 0.333//s, G, =G^=30dB, /^=3GHz, 
/I = 0.1m, crj.=100m^, rQ=290 K, F^=5dB, and Lj = = 1. Also, the 

performance of each P4 code using the corresponding is plotted. For the P4 code 
from Equation (28), / A = 18/3 = 6 , ^c 3 =10- The single P4 

codes have a limited unambiguous range; A^j =6 gives 5^ =3x10* x 10x 6/2 
= 900 m, A ^2 = 8 gives 5,, = 3 x 10* x 10x 8 / 2 = 1200 m, and A^j = 10 gives 5,, = 
3xl0*xl0^®xl0/2= 1500 m. These values are less than achieved using the A = 3 
robust symmetrical residue-P4 phase code sequences, = cM^snsK / 2 = 2150 m. 



Figure 66 . Comparison of the maximum unambiguous range of CW radar system for 
SNRjf^ = 13 dB using robust symmetrical residue-P4 code for m^=3, 

= 4, and =5 with the radar using each P4 code individually for the 
corresponding A = A / A . 
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In conclusion, the robust symmetrical-P4 phase eode sequenee presented 
in this ehapter ean be used to extend the unambiguous range of CW polyphase radar 
systems. The issue eoneerned with the range bin error due to the fluetuation of the 
eompression outputs is also resolved. The integer Gray eode properties ean prevent an 
error in the deteeted robust symmetrieal residues. In the next ehapter, the eoneluding 
results from using the number theoretie transforms from Chapters III to V to extend the 
unambiguous range in CW polyphase radar systems are diseussed. 
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VI. CONCLUDING RESULTS AND RECOMMENDATIONS 


A. EQUATION SUMMARY 

The summary of equations used to eompute the signal and range deteetion 
parameters for eaeh modular phase eode (residue-Frank phase eode, symmetrieal residue- 
P4 phase eode, and robust symmetrieal residue-P4 phase eode) that results in extending 
the unambiguous range are shown in the Table 10. The equations are used for any of the 
eoprime system moduli. 


Table 10. Summary of equations for eaeh modular phase eode. 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Phase sequenee 

m 1 fn j m,k 

, '’■fc ^ 

2) \ 4'" 

<J) 1 ' = -^(RSyy, p'—m] ——m 

^m,k l2m^ ' J 2 

Subeode 

period 

h = cpp / /, 

h = cpp / 

N = cpp/Nf^ 

Signal 

bandwidth (5) 

\/t, 

i/c 


The number of 

subeodes 

) 

^ ^ ^RSNS 


m 

2mN 

Proeessing gain 

(PG) 


m 

2m 

Code period (T) 

m\ 

mt, 

2mNtl 
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Using a set of moduli, we ean extend the maximum deteetion ranges beyond those 
of individual moduli. The equations for ealeulating the dynamie range, the number of 
eode periods within the dynamie range, the range resolution, and the unambiguous range 
are shown in Table 11. 


Table 11. Summary of equations for N - channel modular number system. 


Parameters 

RNS-Frank 

SNS -P4 

RSNS -P4 

Dynamic range 

i^RNS’ ^SNS’^RSNs ) 

N 

i=\ 

\m j N 

mini — Yl mi,+ IT f 

[^1=2 l=j+\ 'J 

(even) 

3 2 15 

— m, H - m, +1 tor 

2 1 2 1 

N = 3 moduli of the 



\\ J \ N ] 
mm -n «;•/+- n nij, 
Pm P 

j> (odd) 

form 2 ^-l, 2 % 2 ^+l 

The number of 

M 

RNS 

9 

^ SNS 


M 

RSNS 

code periods 

m 

m 


2Nm 

(N.J 





Range 





resolution (AR) 

2 

2 


2 

Unambiguous 

cMt^ 

SNS h 


RSNS h 

range (R„) 

2 

2 


2 


In the SNS and RSNS, the dynamic range is not necessarily greater with an 
increasing value of the moduli values. For example, the dynamic range of a 3-channel 
SNS system with moduli [7, 13, 16] is M^ns ~ ^9. On the other hand, the dynamic range 
of a 3-channel RSNS system with moduli [7, 13, 19] is only =55. That is, and 
mj are equal, but the last RSNS modulus is =19, which is larger than the last SNS 
modulus m 3 =16. However, the dynamic range of the RSNS is less than the SNS. 
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B. COMPARISON OF 3-CHANNEL MODULAR NUMBER SYSTEMS WITH 
THE SAME DYNAMIC RANGE 

In Table 12, a comparison of each modular number system having the same 
dynamic range (60), the same subcode period (1 ps), the same unambiguous range (9000 
m), and the same range resolution is shown. Note the bandwidth for the RSNS-P4 system 
is smaller since it is contains N phase values that are equal. 


Table 12. Comparison of each modular phase code in which each has the same 

dynamic range = 60, and t[= \ jus . 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Set; Dynamic range 

60 

60 

60 

Set; Subcode period (ps) 

1 

1 

1 

Bandwidth (MHz) 

1 

1 

0.33 

Unambiguous range (m) 

9,000 

9,000 

9,000 

Range resolution (m) 

150 

150 

150 

Moduli 

[3, 4, 5] 

[7,13,29] 

[3, 5, 7] 

Numbers of phase values 

[3, 4, 5] 

[4, 7, 15] 

[4, 6, 8] 

code period (ps) 

[9,16, 25] 

[7,13,29] 

[18,30, 42] 

Number of code periods 

[6.67, 3.75, 2.4] 

[8.6, 4.6, 2.1] 

[3.3,2, 1.4] 


The RNS compared to the SNS is more efficient in covering since fewer code 
periods are required in the compression. The phase values from the phase detector are 
also less for the same . Comparison with the RSNS shows that it is the most efficient of 

all schemes due to the fewer number of phase values used, the smallest number of code 
periods required in the compression and its Gray code properties that limit the range error 
to one range bin. In terms of complexity, the size of the compressor is longest for the 
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RSNS since the dynamic range incorporates the fewest code period (assuming =1). 
The SNS requires the fewest code periods in the compressor. 

C. COMPARISON OF 3-CHANNEL MODULAR NUMBER SYSTEMS WITH 
THE SAME MODULI 

In Table 13, each number theoretic transform that has the same set of moduli [3 4 
5] with the subeode period equal to 1 ps is shown. 


Table 13. Comparison of each modular number system in which each has the same 

moduli [3 4 5], and = \jus 


Parameters 

RNS-Frank 

SNS-P4 

RSNS-P4 

Set; Moduli 

[3 4 5] 

[3 4 5] 

[3 4 5] 

Dynamie range (units) 

60 

11 

43 

Set: Subcode period (ps) 

1 

1 

1 

Bandwidth (MHz) 

1 

1 

0.33 

Unambiguous range (m) 

9,000 

1,650 

6,450 

Range resolution (m) 

150 

150 

150 

Number of phase values 

[3, 4, 5] 

[2, 3, 3] 

[4, 5, 6] 

Number of code periods 

[6.67, 3.75, 2.4] 

[3.7, 2.8, 2.2] 

[2.4, 1.8, 1.4] 


In the RNS phase code, the extended unambiguous range is maximized due to its 
maximum dynamic range. In the SNS phase code, the number of phase values is lowest 
due to the eonstruction of the SNS integers within the modulus. However, it gives the 
shortest unambiguous range sinee its dynamie range is the smallest. The RSNS gives the 
next largest unambiguous range of the three schemes. The Gray code properties make it 
an attraetive ehoice. 
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D. RANGE DETECTION ERROR CONTROL 


As shown, only the robust symmetrical residue-P4 phase code sequence could 
prevent a range error due to the noise changing the phase of the return signal. This is due 
to the Gray code property that the RSNS has. For the other modular phase codes (residue- 
Frank and symmetrical residue-P4), the range error is large when the target straddles two 
range bins. The higher set of moduli might help control the range error due to its higher 
processing gain. However, increasing the modular number would increase the number of 
phase steps, resulting in increasing receiver’s complexity and a longer code period, 
resulting in higher compression loss and required dwell time on targets. 

E. RECOMMENDATIONS FOR FUTURE WORK 

The modular phase codes used in this thesis were shown to be able to extend the 
unambiguous range in CW polyphase modulation radar systems. The applications of 
these modular phase codes in different types of radar systems would be interesting to 
study. Also, the performance of each moduli set with the different trade-offs should be 
further considered. The Doppler tolerant property that the P4 code has could also be 
investigated for the symmetrical residue-P4 and the robust symmetrical residue-P4 phase 
code sequences. Finally, the hardware implementation has to be examined and a 
simulation that includes hardware imperfections and errors must be conducted. 
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APPENDIX - LIST OF VARIABLES 


There are many variables used throughout this thesis. All of the important 
variables are listed here to provide unambiguous meanings. 


Table 14. List of variables. 


Variable 

Meaning 

a 

Range bin index 

c 

Speed of light (m) 

cpp 

Number of earner eyeles within a subeode 

fc 

Carrier frequeney 

fs 

Sampling frequeney (Hz) 

h 

Integer position in RSNS 

k 

Phase index 

kg 

Boltzman’s eonstant (joule/K) 

m 

Coprime modulus 

m,. 

Coprime moduli 

h 

Subeode period (s) 

B 

Bandwidth (Hz) 


Reeeiver’s input bandwidth (Hz) 

F, 

Reeeiver noise faetor 

G. 

Gain of reeeiving antenna 

G, 

Gain of transmitting antenna 

Frr 

Loss between the radar’s antenna and reeeiver 

Lrt 

Loss between the radar’s transmitter and antenna 

L, 

Two-way atmospherie transmission loss 

M 

Frank frequeney step 

^ RNS 

Dynamie range of the RNS 

^SNS 

Dynamie range of the SNS 

^ RSNS 

Dynamie range of the RSNS 

N 

Number of eoprime modular sequenees 

K 

Number of subeodes within a eode period 

N 

^ ^ cRSNS 

Number of subeodes within a eode period (RSNS) 

Nr 

Number of eode periods returned from target 
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Variable 

Meaning 


Number of eode periods used for modulus m 

N. 

Number of reeeiver referenee eodes 

PG 

Proeessing gain 

PSL 

Peak-to-sidelobe level 

P 

RNS 

Period of a single RNS sequenee 

P 

SNS 

Period of a single SNS sequenee 

P 

-* RSNS 

Period of a single RSNS sequenee 

R 

Range to target (m) 

K 

Unambiguous range (m) 

D 

^u,RNS 

Unambiguous range of RNS (m) 

D 

^u,SNS 

Unambiguous range of SNS (m) 

D 

^u,RSNS 

Unambiguous range of RSNS (m) 

^m.k 

RNS residue kth for modulus m 

^^m,k 

RSNS symmetrieal residue kth for modulus m 

^m,k 

SNS symmetrieal residue kth for modulus m 

SNR,, 

Input signal-to-noise ratio 

SNR,„ 

Output signal-to-noise ratio 

T 

Code period (s) 

T 

nij 

Code period for eoprime modulus m. 

V 

Target veloeity (m/s) 

T 

Roundtrip time delay (s) 

^PW 

Pulse width 

z 

Signal wavelength (m) 

V 

Doppler frequeney (Hz) 

AR 

Range resolution (m) 

CTY 

Target’s radar eross seetion (m ) 

^m,k 

Phase eode for modulus m with index k 
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